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Abstract 

The seemingly simple problem of determining the drag on a body moving through 
a very viscous fluid has, for over 150 years, been a source of theoretical confu- 
sion, mathematical paradoxes, and experimental artifacts, primarily arising from 
the complex boundary layer structure of the flow near the body and at infinity. 
We review the extensive experimental and theoretical literature on this problem, 
with special emphasis on the logical relationship between different approaches. The 
survey begins with the developments of matched asymptotic expansions, and con- 
cludes with a discussion of perturbative renormalization group techniques, adapted 
from quantum field theory to differential equations. The renormalization group 
calculations lead to a new prediction for the drag coefficient, one which can both 
reproduce and surpass the results of matched asymptotics. 
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I. INTRODUCTION TO LOW R FLOW 



A. Overview 



In 1851, shortly after writing down the Navier-Stokes equations, Sir George Gabriel 
Stokes turned his attention to what modern researchers might whimsically refer to as "the 
hydrogen atom" of fluid mechanics: the determination of the d r ag on a sphere or an infinite 
cylinder moving at fixed speed in a highly viscous fluid (jStokesl . Il85ll ). Just as the quantum 
theory of the hydrogen atom entailed enormous mathematical difficulties, ultimately leading 
to the development of quantum field theory, the problem posed by Stokes has turned out 
to be much harder than anyone could reasonably have expected: it took over 100 years 
to obtain a justifiable lowest order approximate solution, and that achievement required 
the invention of a new branch of applied mathematics, matched asymptotic expansions. 
And just as the fine structure of the hydrogen atom's spectral lines eventually required 
renormalization theory to resolve the problems of "infinities" arising in the theory, so too, 
Stokes' problem is plagued by div ergences tha t are, t o a physicist, most naturally resolvec 



by renormalization group theory (IChen et all 1 19961 : iFeynmanl . 



1948 



Gell-Mann and Low 
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1954 



Schwingerl . 1 19481 ; IStuckelberg and Petermannl . Il953l ; iTomonaga 



1948 



Wilson 



1971 



1983|). 

In order to appreciate the fundamental difficulty of such problems, and to expose the 
similarity with familiar problems in quantum electrodynamics, we need to explain how 
perturbation theory is used in fluid dynamics. Every flow that is governed by the Navier- 
Stokes equations only (i.e. the transport of passive scalars, such as temperature, is not 
considered; there are no rotating frames of reference or other complications) is governed by 
a single dimensionless parameter, known as the Reynolds number, which we designate as R. 
The Reynolds number is a dimensionless number made up of a characteristic length scale L, 
a characteristic velocity of the flow U, and the kinematic viscosity v = T]/p, where rj is the 
viscosity and p is the density of the fluid. In the problems at hand, defined precisely below, 
the velocity scale is the input fluid velocity at infinity, Moo, and the length scale is the radius 
a of the body immersed in the fluid. Then the Reynolds number is given by: 

R = ^ (1) 
v 

The Reynolds number is frequently interpreted as the ratio of the inertial to viscous terms 
in the Navier-Stokes equations. For very viscous flows, R — > 0, and so we anticipate that 
a sensible way to proceed is perturbation theory in R about the problems with infinite 
viscosity, i.e. R = 0. In this respect, the unwary reader might regard this as an example 
very similar to quantum electrodynamics, where the small parameter is the fine structure 
constant. However, as we will see in detail below, there is a qualitative difference between 
a flow with R = and a flow with R — > 0. The fundamental reason is that by virtue of 
the circular or spherical geometry, the ratio of inertial to viscous forces in the Navier-Stokes 
equations is not a constant everywhere in space: it varies as a function of radial distance r 
from the body, scaling as O (Rr/a). Thus, when R = 0, this term is everywhere zero; but 
for any non-zero R, as r/a — > oo the ratio of inertial to viscous forces becomes arbitrarily 
large. Thus, inertial forces can not legitimately be regarded as negligible with respect to 
viscous forces everywhere: the basic premise of perturbation theory is not valid. 

Perturbation theory has to somehow express, or manifest, this fact, and it registers its 
objection by generating divergent terms in its expansion. These divergences are not physical, 
but are the perturbation theory's way of indicating that the zeroth order solution — the point 
about which perturbation theory proceeds — is not a correct starting point. The reader 
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might wonder if the precise nature of the breakdown of perturbation theory, signified by 
the divergences, can be used to deduce what starting point would be a valid one. The 
answer is yes: this procedure is known as the perturbative renormalization group (RG), 
and we will devote a significant fractio n of this article to expounding this strategy As 



most readers will know, reno rmalization (IFevnman 



and renormali za tion g roup (jGell-Mann and Lowl . 



Wilsonl . 



1971a 



1948 



1954 



Schwingerl. 



194 



Tomonaea . 



Stuckelberg and Petermann 



1948) 



1953 



19831 ) techniques in quantum field theories have been stunningly successful. 



In the most well-controlled case, that of quantum electrodynamics, the smallness of the 
fine structure constant allows agreem ent of perturbative ca lculations with high-precision 
measurements to 12 significant figures (IGabrielse et all 120061 ). Do corresponding techniques 
work as well in low Reynolds fluid dynamics, where one wishes to calculate and measure 
the drag Cd (defined precisely below)? Note that in this case, it is the functional form in 
R for the drag that is of interest, rather than the drag at one particular value of R, so the 
measure of success is rather more involved. Nevertheless, we will see that calculations can 
be compared with experiments, but there too will require careful interpretation. 

Historically a different strategy was followed, leading to a set of techniques known generi- 
cally as singular perturbation theory, in particular encompassing boundary layer theory and 
the method of matched asymptotic expansions. We will explain these techniques, developed 
by mathematicians starting in the 1950's, and show their connection with renormalization 
group methods. 

Although the calculational techniques of matched asymptotic expansions are widely 
regarded as representing a systematically firm footing, their best results apply only to 
infinitesimally small Reynolds number. As shown in Figure [TJ the large deviations be- 
tween theory and experiment for R ~ 0.5 demonstrate the need for theoretical predictions 
which are more robust for small but non-infinitesimal Reynolds numbers. Ian Proudman, 
who, in a tour de force helpe d obtain the first matched asymptotics result for a sphere 



(jProudman and Pearson! 



19571 ). expressed it this way: "It is therefo re particularly disap 



point ing that the numerical 'convergence' of the expansion is so poor." (IChester and Breach! . 



19691 ) In spite of its failings, Proudman's solution from 1959 was the first mathematically 



rigorous one for flow past a sphere; all preceding theoretical efforts were worse. 

Further complicating matters, the literature surrounding these problems is rife with 
"paradoxes", revisions, ad-hoc justifications, disagreements over attribution, mysterious fac- 
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1.5 



1.5 



-Eqn. [63] 
Tritton 
Jayaweera 




□ 



• Dennis 
X Maxworthy 
+ Le Clair 
-Eqn. [55] 




FIG. 1 ( Color online) Comparing experiment with "state of the art" theor etical predictions for 



1971 



a sphere dJavaweera and Mason 



Le Clair and Hamielec 



1965 



1970|; 



Tritton 



Maxworthy 



1959) (right) and a cylinder ([Dennis and Walkerl . 



196.51 ) (left). 



tors of two, conflicting terminology, non-standard d efinitions, and la nguage barriers. Even a 
recent article attempting to resolve this quagmire (ILindgrenl . Il999l ) contains an inaccuracy 
regarding publication dates and scientific priority. This tortured history has left a wake of 
experiments and numerical calculations which are of widely varying quality, although they 
can appear to agree when not examined closely. For example, it turns out that the finite size 
of experimental systems has a dramatic effect on measurements and simulations, a problem 
not appreciated by early workers. 

Although in principle the matched asymptotics results can be systematically extended 
by working to higher order, this is not practical. The complexity of the governing equations 
prohibits further improvement. We will show here that techniques based on the renormal- 
ization group ameliorate some of the technical difficulties, and result in a more accurate drag 
coefficient at small but non-infinitesimal Reynolds numbers. Given the historical importance 
of the techniques developed to solve these problems, we hope that our solutions will be of 
general methodological interest. 

We anticipate that some of our readers will be fluid dynamicists interested in assessing 
the potential value of renormalization group techniques. We hope that this community 
will see that our use of the renormalization group is quite distinct from applications to 
stochastic problems, such as turbulence, and can serve a different purpose. The second 
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group of readers may be physicists with a field theoretic background, encountering fluids 



problems for the first time, perhaps in unconventional settings, such as heavy ion co 



and QCD (lAckermann et al. 



Hirano and Gyulassyl . 



2001 



Baier et al. 



2006 



20061 ) or 2D electron gases (jEavesl . 



Csernai et al 



1998 : 



2006, 



Stone 



2005 



Heniz 



lisions 



2005 



1990). We hope that 



this review will expose them to the mathematical richness of even the simplest flow settings, 
and introduce a familiar conceptual tool in a non-traditional context. 

This review has two main purposes. The first purpose of the present article is to attempt 
a review and synthesis of the literature, sufficiently detailed that the subtle differences 
between different approaches are exposed, and can be evaluated by the reader. This is 
especially important, because this is one of those problems so detested by students, in which 
there are a myriad of ways to achieve the right answer for the wrong reasons. This article 
highlights all of these. 

A second purpose of this article is to review the use of renormalization group techniques 
in the context of singular perturbation theory, as applied to low Reynolds number flows. 
These techniques generate a non-trivial estimate for the functional form of Cd(R) that 
can be sensibly used at moderate values of R ~ 0(1), not just infinitesimal values of R. 
As R — > 0, these new results reduce to those previously obtained by matched asymptotic 
expansions, in particular accounting for the nature of the mathematical singularities that 
must be assumed to be present for the asymptotic matching procedure to work. 

Renormalization group techniques were originally developed in the 1950's to extend and 
improve the perturbation theory for quantum electrodynamics. During the late 1960's 



and 1970's, renormali zation group te c hniqu e s famou s 



of phase transitions (IL. P. Kadanofi . 



1966 



Widoml . 



y found appli c ation in the problem 



1963 



Wilsonl . Il971al ). During the 



1990's, renormalization group techniques were developed for ordinary and partial differ- 
ential equations, at first for the analysis of no nequilibrium (but deterministic) pr o blems 



which exhibited anomalous scaling exponents (IChen et al. 



1991 



Goldenfeld et al 



and subsequently 



1995 



Chen et al. 



or the related problem of travelling wave selection (jChen and Goldenfeldl . 



1990) 



1994a 



c). The most recent significant development of the renormalization 



group — an d the one that concerns us here — was t he application to singular perturbation 



problems (IChen et al. 



1994b 



19961 ). The scope of (IChen et all Il996l ) encompasses bound- 



ary layer theory, matched asymptotic expansions, multiple scales analysis, WKB theory, and 
reductive perturbation theory for spatially-extended dynamical systems. We do not review 
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all these developments here, but focus only on the issues arising in the highly pathological 
singularities characteristic of low Reynolds number flows. For a pedagogical i ntroduction to 
renormalization group techniques, we refer the reader to (jGoldenfeldl . Il992l ). in particular 
Chapter 10 which explains the connection between anomalous dimensions in field theory 



and similarity solutions of partial differential equations. We mention a 
niques discussed here haye also been the subject of r i gorous analysis ( 



Bricmont and Kupiainen 



Moise and 



2002 



Ziane, 



emam 



2000 



995 



Bricmont et al. 



Moise and Ziane 



1994 



2001 



Lan and 



Petcu et al 



AW 



so that the RG tech- 



Blomker et al. 



200 



2005 



Moise et al 



2002; 



1998; 



Wirosoetisno et al. 



in cavitation ( Josser and, 



Iguchi et al. 



2000T) in ot h er co ntexts of fluid dynamics, and have also found app 
1999) and cosmological fluid dynam i cs (|Belinchon et al. 



1998 : 



Nambul . 



2000 



2002 



ication 



Nambu and Yamaguchil . 



2002 ; 



1999|). 



This review is organized as follows. After precisely posing the mathematical problem, we 
review all prior theoretical and experimental results. We identify the five calculations and 
measurements which are accurate enough, and which extend to sufficiently small Reynolds 
number, to be useful for evaluating theoretical predictions. Furthermore, we review the his- 
tory of all theoretical contributions, and clearly present the methodologies and approxima- 
tions behind previous solutions. In doing so, we eliminate prior confusion over chronology 
and attribution. We conclude by comparing the best experimental results with our new, 
RG-based, theoretical prediction. This exercise makes the shortcomings that Proudman 
lamented clear. 



B. Mathematical formulation 

The goal of these calculations is to determine the drag force exerted on a sphere and on 
an infinite cylinder by steady, incompressible, viscous flows. The actual physical problem 
concerns a body moving at constant velocity in an infinite fluid, where the fluid is at rest 
in the laboratory frame. In practice, it is more convenient to analyze the problem using an 
inertial frame moving with the fixed body, an approach which is entirely equivalent. 1 

Flow past a sphere or circle is shown schematically in Figure [2j The body has a char- 
acteristic length scale, which we have chosen to be the radius (a), and it is immersed in 



1 Near ly all workers, beginning with Stokes (jStoked . 1 185 If ) . use this approach, which Lindgren (jLindgren . 



1999) refers to as the "steady" flow problem. 
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1 \ 




r 











FIG. 2 (Color online) Schematic for flow past a sphere or cylinder. 



Quantity 


Description 


r 


Coordinate Vector 


u{r) 


Velocity Field 


P 


Fluid Density 


p{r) 


Pressure 


V 


Kinematic Viscosity 


a 


Characteristic Length of Fixed Body 




The Uniform Stream Velocity 



TABLE I Quantities needed to characterize low R flow past a rigid body. 



uniform stream of fluid. At large distances, the undisturbed fluid moves with velocity Mqo. 

The quantities shown in Tabled characterize the problem. We assume incompressible flow, 
so p = const. The continuity equation (Eqn. [2]) and the time-independent Navier-Stokes 
equations (Eqn. [3]) govern steady-state, incompressible flow. 

V-m = (2) 

(u-Vn) = -— + z/V 2 m (3) 
P 

These equations must be solved subject to two boundary conditions, given in Eqn. |U First, 
the no-slip conditions are imposed on the surface of the fixed body (Eqn. Haj). Secondly, the 
flow must be a uniform stream far from the body (Eqn. Ebl) . To calculate the pressure, one 
also needs to specify an appropriate boundary condition (Eqn. He]) , although as a matter 
of practice this is immaterial, as only pressure differences matter when calculating the drag 
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Dimensionless Quantity 


Definition 




f/a 


u*(f) 


«(r)/|«oo| 


p*(f) 


a p(f)/p 1/ Uoo 


V* 


a V 



TABLE II Dimensionless variables. 



coefficient. 

u{r) = r G {Surface of Fixed Body} (4a) 

lim u(r) = (4b) 

lim p(r) = (4c) 

It is convenient to analyze the problem using non-dimensional quantities, which are de- 
fined in Table [Til When using dimensionless variables, the governing equations assume the 
forms given in Eqns. [5]and[6l where we have introduced the Reynolds Number, R = [u^a/v, 
and denoted scaled quantities by an asterisk. 

V* • u* = (5) 

R(u* ■ V*)u* = -W + V*V (6) 

The boundary conditions also transform, and will later be given separately for both the 
sphere and the cylinder (Eqns. [TH [TUl) . Henceforth, the * will be omitted from our nota- 
tion, except when dimensional quantities are explicitly introduced. It is useful to eliminate 
pressure from Eqn. [6] by taking the curl and using the identity V x Vp = 0, leading to 

(u ■ V)(V x u) - ((V x u) ■ u) = ^V 2 (V x u) (7) 

R 

1. Flow past a cylinder 

For the problem of the infinite cylinder, it is natural to use cylindrical coordinates, r = 
(r, 6, z). We examine the problem where the uniform flow is in the x direction (see Figure 
E]). We will look for 2-d solutions, which satisfy d z u = 0. 
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Since the problem is two dimensional, one may reduce the set of governing equations 
(Eqns. [5] and |6| to a single equation involving a scalar quantity, the Lagrangian stream 
function, usually denoted ip(r,9). It is defined by Eqn. [SI 2 



u. 



I dip 
r 09 



Uq = - 



dip 
dr 







(8) 



This definition guarantees that equation (J5]) will be satisfied (IGoldsteinl . Il929l ) . Substi- 
tuting the stream function into equation (GO), one obtains the governing e quation (Eqn . 



9]). Here we follow the co mpact notation of Proudman and Pearson (iHinchl . 



Proudman and Pearson 



1991 



1953). 



Rd{iP,V 



where 



d 2 



+ 



r d(r, 9) 



Id Id 2 
+ 



(9) 



dr 2 r dr r 2 d9 2 

The boundary conditions which fix u(r) (Eqns. HaJ Hbl) also determine ip(r,9) up to an 
irrelevant additive constant. 3 Eqn. [10] gives the boundary conditions expressed in terms of 
stream functions. 



ip(r = 1 
dip(r, 9) 



dr 



lim 



iP(r,9) 






sin(0) 



(10a) 
(10b) 

(10c) 



r— >oo T 

To calculate the drag for on a cylinder, we must first solve Equation [9] subject to the 
boundary conditions given in Eqn. [101 



2. Flow past a sphere 

To study flow past a sphere, we use spherical coordinates: r = (r,9,(p). We take the 
uniform flow to be in the z direction. Consequently, we are interested in solutions which are 
independent of <p, because there can be no circulation about the z axis. 



2 



Although many authors pref er to solve the vector equations, we follow Proudman and Pearson 



( Proudman and Pearson! 119571 ) 



3 The constant is irrelevant because it vanishes when the derivatives are taken in Eqn. [H 
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Since the problem has axial symmetry, one can use the Stokes ' stream function (or Stokes' 
current function) to reduce Eqns. [5] and [6] to a single equation. This stream function is 
defined through the following relations: 

1 1 



V<j> 







(11) 



r 2 sin 9 r sin 9 

These definitions guarantee that Eqn. [5] will be satisfied. Substituting Eqn. [TTl into Eqn. [3, 
one obtains the governing equation for ip(r, 9) f Proudman and Pearsorl 1957 ): 



- a(r, ji) r z 



(12) 



In this equation, 



/i = cos 9 

c?_ 1 

dr 2 



D 2 = 
L = 



+ 



fi 2 d 2 



d 1 d 
+ 



1 — fi 2 dr r d\i 



Here we follow the notation of Proudman and Pearson (IProudman and Pearson 



Other authors, such as Van Dyke (IVan Dykd . Il975l ) and Hinch (IHinchl . Il99ll ). write their 



19571 ). 



stream function equations in an equivalent, albeit less compact, notation. 

As in the case of the cylinder, the boundary conditions which fix u[f) (Eqns. HaJ l4b|) 
determine ip up to an irrelevant additive constant. The transformed boundary conditions 
are given by Eqn. \l£ 



if>(r = l,fi) 
dip(r, /i) 



dr 



lim 

r— »oo 



r=l 






\ (I-. 2 ) 



(14a) 
(14b) 

(14c) 



In this paper, we obtain approximate solutions for Eqn. [9] (subject to Eqn. [10]), and Eqn. 
(subject to Eqn. [H]). These solutions are then used to calculate drag coefficients, which 
we compare to experimental results. 



3. Calculating the drag coefficient 

Once the Navier-Stokes equations have been solved, and the stream function is known, 
calculating the drag coefficient, Co, is a mechanical procedure. We follow the methodology 

12 



with the work done by Kaplun (IKaplun 



described by Chester and Breach ([Chester and Breach! . Il969f ) . T his analysis is consistent 



19571 ) and Proudman (jProudman and Pearson 



19571 ). although these authors do not detail their calculations. 



This methodo 



Tomotika (IQseen 



qgy is significantly different from that employed by other workers, such as 



1910 



Tomotika and Aoi 



19501 ) . Tomotika calculates Cp approximately, 



based on a linearized calculation of pressure. Although these approximations are consistent 
with the approximations inherent in their solution of the Navier-Stokes equations, they are 
inadequate for the purposes of obtaining a systematic approximation to any desired order 
of accuracy. 

Calculating the drag on the body begins by determining the force exerted on the body 
by the moving fluid. Usin g dimensional variables, the force per unit area is given by 
([Landau and Lifschita . Il999l ): 

Pi = -&ikn k (15) 
Here is the stress tensor, and n is a unit v ector normal to the surface. For an incom- 



pressible fluid, the stress tensor takes the form (ILandau and Lifschitz 



1999|): 



dvi dv k \ 
dxij 



dx k 



(16) 



7] is the dynamic viscosity, related to the kinematic viscosity by rj = up. The total force is 
found by integrating Eqn. [15] over the surface of the solid body. We now use these relations 
to derive explicit formula, expressed in terms of stream functions, for both the sphere and 
the cylinder. 

a. Cylinder In the case of the cylinder, the components of the velocity field are given 
through the definition of the Lagrangian stream function (Eqn. [8]). Symmetry requires that 
the net force on the cylinder must be in the same direction as the uniform stream. Because 
the uniform stream is in the x direction, if follows from Eqns. [TBI and [TBI that the force 4 on 
the cylinder per unit length is given by: 



Fx = j)(o~rr cos 9 — o~ r Q sin 9)ds 



(17) 



2tt 



. cos 9 



2n 



-p + 2?7 



a r e sin 

dv. 



r d9 



dr 



r 39 dr 



r 



sm9 )r d9 



The form of cr^ in cylindrical coordinates is given is Landau (jLandau and Lifschita . Il999f ) . 



13 



The drag co efficient for an inf i nite c y 



authors (e.g., (ILagerstrom et al 



1967 



inder is defin ed as Cd = F^et/ p\uoo\ 2 a. Note that 



Trittonl . 



19591 )) who define the Reynolds number 



based on diameter nonetheless use the same definition of Cd, which is based on the radius. 
For this problem, F^ et = F%, as given by Eqn. [TT1 Introducing the dimensionless variables 
defined in Table [TT] into Eqn. [T71 we obtain Eqn. [TSJ Combining this with the definition of 
Cd, we obtain Eqn. 



p\u c 



R 



2n 



-p(r, 9) + 2 



8u r 



Or 



n .1 du r duf, 
cos 9— I — tztt + 



r 89 dr 



sin 9 I r d9 



Cd ~r 



2tt 



I JO 



-p(r, 9) + 2 



dr 



ldu r , dUf, 

COS 9 — I -^ttt + 



r 9^ dr 



Ug 

r 



— \ sin9 \ r d9 



r=l 



;is) 



(19) 



J r=l 



To evaluate this expression, we must first derive p(r, 9) from the stream function. 
The pressure can be determined to within an irrelevant addi tive constant by in t egrat 
ing the 9 component of the Navier-Stokes equations (Eqn. [6]) ( Chester and Breach . 



Landau and Lifschit2. 



1969 ; 



19991 ). The constant is irrelevant because, in Eqn. [191 C cos9d9 



0. Note that all gradient terms involving z vanish by construction. 



p(r, 9) 



J -R([u-V)u e + 



u r ue 



f V 2 + ^ ^ Ur 

r 2 09 



Ug 



d9 



(20) 



Given a solution for the stream function ip, the set of dimensionless Eqns. [U [191, and 
uniquely determine Cd for a cylinder. However, because the velocity field satisfies no-slip 
boundary conditions, these general formula often simplify considerably. 

For instance, consider the class of stream functions which meets the boundary conditions 
(Eqn. [TO]) and can be expressed as a Fourier sine series: ip(r, 9) = J^^Li /«( r ) sinn#. Using 
the boundary conditions it can be shown that, for these stream functions, Eqn. [T9l reduces 
to the simple expression given by Eqn. 



C 



71 



D 



Mr) 



(21) 



r=l 



R Vdr 3 ' 
b. Sphere 

The procedure for calculating Cd in the case of the sphere is nearly identical to that for 
the cylinder. The components of the velocity field are given through the definition of the 
Stokes' stream function (Eqn. [TT]) . As before, symmetry requires that any net force on the 
cylinder must be in the direction of the uniform stream, in this case the z direction. 
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From Eqn. [151 the net force on the sphere is given by Eqn. 



[a rr cos 9 — a r e sin 9)ds 



(22) 



2tt 



(a rr cos 9 — o t q sin I 



r 2 sin 9 d9 



For the sphere, the drag coefficient is defined as Cd = F^ et /p\u c 



12^2 



Often the drag 



coefficient is given in terms of the Stokes' Drag, D$ = 6irp\u 00 \ai' = 67rp|w 00 | 2 a 2 /.R. In these 
term s, Cd = F^ et 6ir/ DsR- If F^ et = Ds, Cd = Qir/R, the famous result of Stokes (jstokes . 
18511) . 



Not all authors fol 



([Goldstein 



1929 



ow Stokes' original definition of Cr>- For instance, S . Gold stein 



19651 ) and H. Liebster ( [Liebster 



1927 



Liebster and Schiller 



1924 ) de 



fine Cd using a factor based on cross-sectional C 



Goldstein 
D 



Cd2/tt. These authors 



also define R using the diameter of the sphere rather than the radius. S . Dennis, defines Cn 



simi 



197lh . 



arly to Goldstein, but without the factor of two: C^ 



Dennis 



Co/it (IDennis and Walkerl . 



Using the form of Eqn. [TBI given in Landau (ILandau and Lifschita . ll999l ) and introducing 
the dimensionless variables defined in Table [III into Eqn. [221 we obtain Eqn. [231 Combining 
this with the definition of Cd-, we obtain Eqn. [2] 



D s 

F~ z = — 
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(23) 



(24) 



As with the cylinder, the pressure can be determined to within an irrelevant addi- 
tive constant by int e grating the 9 component o f the Navier-Stokes equations (Eqn. EJ) 



( [Chester and Breach! , 
ing (J) must vanish. 

P(r,9) = 



1969 



Landau and Lifschita . Il999l ). Note that gradient terms involv- 



J -R^u-V)u e + ^-)+V 2 u e + 



2 du r 



U0 



d9 



(25) 



r 2 09 r 2 sin 9 2 

Given a solution for the stream function if), the set of dimensionless Eqns. HH [231 an d 
uniquely determine Cd for a sphere. 

As with the cylinder, the imposition of no-slip boundary conditions considerably simpli 
fies these general formula. In particular, consider stream functions of the form ip(r,9) = 
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Y^n=i fn(r)Q n (cos9), where Q n (x) is defined as in Eqn. [46] If these stream functions satisfy 
the boundary conditions, the drag is given by Eqn. [26] 

CD = jk(- 2f,;{r) + f ' ; ' {r) )r^ (26) 

c. A subtle point 

When applicable, Eqns. [21] and [26] are the most convenient way to calculate the drag 
given a stream function. They simply require differentiation of a single angular term's 
radial coefficient. However, they only apply to functions that can be expressed as a series of 
harmonic functions. Moreover, for these simple formula to apply, the series expansions must 
meet the boundary conditions exactly. This requirement implies that each of the functions 
fi(r) independently meets the boundary conditions. 

The goal of our work is to derive and understand approximate solutions to the Navier- 
Stokes' equations. These approximate solutions generally will not satisfy the boundary 
conditions exactly. What — if any — applicability do Eqns. [21] and [26] have if the stream 
function does not exactly meet the boundary conditions? 

In some rare cases, the stream function of interest can be expressed in a convenient 
closed form. In these cases, it is natural to calculate the drag coefficient using the full set 
of equations. However we will see that the solution to these problems is generally only 
expressible as a series in harmonic functions. In these cases, it actually preferable to use the 
simplified equations [21] and [26] 

First, these equations reflect the essential symmetry of the problem, the symmetry im- 
posed by the uniform flow. Eqns. [21] and [26] explicitly demonstrate that, given an exact 
solution, only the lowest harmonic will matter: Only terms which have the same angular 
dependence as the uniform stream will contribute to the drag. By utilizing the simplified 
formula for Co as opposed to the general procedure, we effectively discard contributions 
from higher harmonics. This is exactly what we want, since these contributions are artifacts 
of our approximations, and would not be present in an exact solution. 

The contributions from inaccuracies in how the lowest harmonic meets the boundary 
conditions are more subtle. As long as the boundary conditions are satisfied to the accuracy 
of the overall approximation, it does not matter whether one uses the full-blown or simpli- 
fied drag formula. The drag coefficients will agree to within the accuracy of the original 
approximation. 
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In general, we will use the simpli fied formula. This is the approach taken explicitly by 



many matched asympt otics workers ( 



Chester and Breach!. 



1957 



1969 



Van Dyke 



Skinner 



19751 ). and implic- 



19751 ). It should be noted 



itly by other workers (jProudman and Pearson! 
that these workers only use the portion 5 of their solutions which can exactly meet the as- 
sumptions of the simplified drag formula. However, as we will subsequently discuss, this is 
an oversimplification. 



II. HISTORY OF LOW R FLOW STUDIES 
A. Experiments and numerical calculations 

Theoretical attempts to determine the drag by solving the Navier-Stokes' equations have 
been paralleled by an equally intricate set of experiments. In the case of the sphere, ex- 
periments usually measured the terminal velocity of small falling spheres in a homogeneous 
fluid. In the case of the cylinder, workers measured the force exerted on thin wires or fibers 
immersed in a uniformly flowing viscous fluid. 

These experiments, while simple in concept, were difficult undertakings. The regime 
of interest necessitates some combination of small objects, slow motion, and viscous fluid. 
Precise measurements are not easy, and neither is insuring that the experiment actually 
examines the same quantities that the theory predicts. All theoretical drag coefficients 
concern objects in an infinite fluid, which asymptotically tends to a uniform stream. Any 
real drag coefficient measurements must take care to avoid affects due to the finite size of 
the experiment. Due to the wide variety of reported results in the literature, we found it 
necessary to make a complete survey, as presented in this section. 



1. Measuring the drag on a sphere 



As mentioned, experiments measuring the drag on a sphere at low Reynolds number were 
intertwined with theoretical developments. Early experiments, which essentially c onfirmed 



Stokes' law as a rea sonable ap proximation, inc lude those of Allen 



( I Arnold H. D.I . Il91ll ) , Williams (IWiUiamsl . Il915l ) , and Wieselsberger (jWieselsbergerl . 119221 ) 



Allen, 



19001) . Arno 



To be precise, they use only the Stokes' expansion, rather than a uniform expansion. 
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FIG. 3 (Color online) Early measurements of the drag on a sphere ([Goldstein 



1963 ). 



The next round of experimen ts were done in the 1920s, motivated by the theoretical 
advances b e gun b y C. W. Oseen ioseen . 1910 ). These experimentalists inc luded Schmeidel 



(ISchmiedell . Il928l ) and Liebster ( ILiebsterl . 



1927 



Liebster and Schiller 



192J). T he results of 



Allen , Liebster, and Arnold were analyzed, collated, and averaged by Castleman (jCastlemanl . 



19251 ). whose paper is often cited as a summary of prior experiment s. The sta t e of a ffairs 



after t his work is w ell summarized in plots given by Goldstein (p. 16) ([Goldstein! Il965l ) . and 
Perry (jPerryl . Il950l ). Figure [3] shows Goldstein's plot, digitized and re-expressed in terms of 
the conventional definitions of Cp and R. 

Figure [3] shows the experimental data at this point, prior to the next theoretical devel- 
opment, matched asymptotics. Although the experimental data seem to paint a consistent 
portrait of the function Cd(R), in reality they are not good enough to discriminate between 
different theoretical predictions. 

Finite geometries cause th e most significant exp erimental errors for these measurements 



( jLindgrenl . 



1999 



Maxworthyl . 



1965 



Trittonl . Il988l ). Tritton notes that "the container di- 



ameter must be more than one hundred times the sphere diameter for the error to be less 
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than 2 percent", and Lindgren estimates that a ratio of 50 between the container and sphere 
diameters will result in a 4% change in drag force. 

In 1961, Fidleris et al. ex perimentally studied th e effec ts of finite container size on 



drag coefficient measurements (IFidleris and Whitmord . 



196ll ). They concluded that there 



were significant finite size effects in previous experiments, but also proposed corrections 



to compensat e for ear 



experiments (jLindgrenl . 



i er exp erimental limitations. Lindgren also conducted some related 
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T. Maxworthy also realized this problem, and undertook experiments which could be 
used to evaluate the more precise predictions of matched asymptotics theories. In his own 
words, 

From the data plotted in Goldstein or Perry, it would appear that the presently 
available data is sufficient to accurately answer any reasonable question. How- 
ever, when the data is plotted 'correctly'; that is, the drag is non-dimensionalized 
with respect to the Stokes drag, startling inaccuracies appear. It is in fact im- 
possible to be sure of the drag to better than ±20% ... The difficulties faced 
by previous investigators se emed to be mainly due to an inability to accurately 
compensate for wall effects (iMaxworthyl . Il965l ) . 



Maxworthy refined the falling sphere technique to produce the best experimental mea- 
surements yet — 2% error. He also proposed a new way of plotting the data, which removes 
the divergence in Eqn. [2H (as R—*0). His approach makes clear the failings of earlier 
measurements, as can be seen in Figure HI where the drag measurements are normalized by 
the Stokes drag, Cf okes = 6tt/R. 

In Maxworthy's apparatus, the container diameter is over 700 times the sphere diameter, 
and does not contribute significantly to experimental error, which he estimates at better 
than 2 percent. Note that the data in Figure H] are digitized from his paper, as raw data are 
not available. 

This problem also attracted the attention of atmospheric scientists, who real- 
ized its significance in cloud physics, where " cl oud drops may well be approxi- 



(|Beard and Pruppacher 



mated by rigid spheres." (IPruppacher and Le Clairl. Il970l) In a series of papers (e.g. 



1969 



Pruppacher and Steinberger 



Le Clair and Hamielec 



197(1 



Pruppacher and Le Clair 



1970 ; 



19681 )). H.R. Pruppacher and others undertook numerical and 
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FIG. 4 (Color online) Maxworthy's accurate measurements o f the drag on a sphere (jMaxworthvl . 



19651 ) contrasted with previous experiments (jGoldsteinl . 



1963 ). 



experimental studies of the drag on the sphere. They were motivated by many of the same 
reasons as Maxworthy, because his experiments covered only Reynolds numbers between 
0.4 and 11, and because "Maxworthy's experimental setup and procedure left considerable 
room for improvement" (IPruppacher and Steinbergerl . Il968l ). 

Their results included over 220 measurements, which they binned and averaged. They 
presented their results in the form of a set of linear fits. Adopting Maxworthy's normaliza- 
tion, we collate and summarize their findings in Eqn. [271 



c R 

07T 



0.102 (2R) u - yb0 0.005 < R < 1.0 
0.115 (2i?) - 802 1.0 < R< 20 



(27) 



0.189 (2R) om2 20 < R < 200 



Unfortunately, one of their later papers includes the following footnote (in our notation): 
"At R < 1 the most recent values of CdR/6tt — 1 (Pruppacher, 19 69, unpublished) tended 
to be somewhat higher than those of Pruppacher and Steinberger." (ILe Clair and Hamieled . 



19701 ) Their subsequent papers plot these unpublished data as "experimental scatter." As 
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jetter agreement with both M axworthy's measurements 



1970J), it makes us question the 



the unpublished data are in much 

and their own numerical analysis (ILe Clair and Hamieled . 
accuracy of the results given in Eqn. [271 

There are many ot her numerical calcu l ations of the drag c oefficient for a sphere 
inclu di ng: Dennis 



970: 



Dennis and Walker 



Pruppacher and Le Clair 



( IRimon and Cheng 



Le Clair ( 


Le Clair and Hamielec, 


(Hamielec et al. 




1967 


), 1 


Pinion 


and Kawaguti 


Kawaguti. 


1950). 



Most of these results are not useful either because of large errors (e.g., Jenson), or because 
they study ranges of Reynolds number which do not include R < 1. Many numerical studies 
examine only a few (or even just a single) Reynolds numbers. For the p urposes of comparing 



theor etical predictio ns of Cp at low Reyno 
19711 ) and Le Clair (ILe Clair and Hamieled . 



ds nu mber, only Dennis (IDennis and Walker! . 



19701 ) have useful calculations. Both of these 



papers report tabulated results which are in very good agreement with both each other and 
Maxworthy; at R = 0.5, the three sets of results agree to within 1% in Co, and to within 
10% in the transformed variable, CdR/Qk — 1- The agreement is even better for R < 0.5. 

Figure [5] shows all relevant experimental and numerical results for the drag on a 
sphere. Note the clear disagreement between Pruppacher's results (Eqn. [27]) . and all 
of the other results for R < 1 — including Le Clair and Pruppacher's numerical results 



(ILe Clair and Hamieled . Il970l ). This can be clearly seen in the inset graph. Although 
Pruppacher's experiment results do agree very well with other data for larger values of R 
(R > 20), we will disregard them for the purposes of evaluating theoretical predictions at 
low Reynolds number. 

It should also be noted that there is a community of researchers interested in sedimenta- 
tion and settling velocities who have studied the drag on a sphere. In a contribution to this 
l iterature, Brown reviews all of the authors discussed here, as he tabulates Cr> for R < 5000 
( jBrown and Lawlerl . 120031 ). His report addresses a larger range of Reynolds numbers and 
he summarizes a nu mber of experiments not trea ted here. His methodology is to apply the 
Fidleris' correction (IFidleris and WhitmoreL Il96lh to previous experiments where tabulated 
experimental data was published. 6 While this yields a reasonably well-behaved drag coeffi- 
cient for a wide range of Reynolds numbers, it is not particularly useful for our purposes, as 



Brown incorrectly reports Dennis' work (|Dennis and Walkerl . Il97ll ) as experimental 
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ID 




? IG. 5 (Color online) A 



( Dennis and Walkerl . 



1971 



summary of experimental and numerica 



Le Clair and Hamielec 



1970; 



Maxworthv 



studi es of Cd for a sphere 



19651 ). 



less accurate work obfuscates the results of the most precise experiments near R = 0. It also 
does not include numerical work or important results which are only available graphically 
(e.g., Maxworthy (iMaxworthyl . Il965l )). 



2. Measuring the drag on a cylinder 



Experiments designed to measure the drag on an infinite cylinder in a uniform fluid came 
later than those for spheres. In addition to being a more difficult experiment — theoretical 
calculations assume the cylin der is infinite — there were no theoretical predictions to test 



before Lamb's result in 1911 ( Lamb H. 



i9iih . 



In 1914, E. F. Relf conducted the first experiments (jRelfl , I1914J ). These looked at the 
force exerted on long wires in a fluid. Relf measured the drag down to a Reynolds number 
of about ten. In 1921, Wieselberger measured the drag at still lower Reynods number, 



reachin g R = 2.11 by 



stream (jWieselsbergerl . 



pokin g at the deflection of a weight suspended on a wire in an air 



192l|). 
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These experiments, combined with others ([Goldstein! . Il965l ; iLinkd . Il93ll ) at higher 
Reynolds number, characterize the drag over a range of Reynolds numbers (see Goldstein, 
pg. 15). However, they do not probe truly small Reynolds numbers (R <C 1), and are of 
little use for evaluating theories which are only valid in that range. Curiously, there are 
no shortage of claims otherwise, such as Lamb, who says "The formula is stated to be in 
good agreem ent with experiment for sufficiently small values of U oc a/i'; see Wieselsberger" 



( ILambl . 



1932|). 



In 1933, Thorn measured the "pressure drag", extending observations down to R = 1.75. 
Thorn also notes that this Reynolds number is still too high to compare with calculations: 
"Actually, Lamb's solution only applies to values of R less than those shown, in fact to values 
much less th an unity, but evidently in most cases the experimental results are converging 
with them." ( Ttioml L933I ) 

In 1946, W hite undertook a series of measurements, which were flawed due to wall effects 



( 1 Whitd . Il946l ). The first high quality expe r iment s which measured the drag at low Reynolds 



number were done by R. K. Finn (Finn. 1 953^. His results, available only in graphical 



form, are reproduced in Figure [7J While vastly superior to any previous results, there is 
considerable scatter in Finn's measurements, and they have largely been surpassed by later 
experiments. 

Tritton, in 1959, conducted experiments which re ached a Reyno lds number of R = 0.2, 



and also filled in some gaps in the R — Co curve (iTrittonl . 119591 ). Tritton estimates his 
accuracy at ±6%, and compares his results favorably to previous work, commenting that, 
"Probably the lowest R points of the other workers were stretching their techniques a little 
beyond their limits." Tritton is also the first author to give a discussion of systematic errors. 7 
Tritton's results are shown in Figure [6j All of his data are available in tabular form. 

Maxworthy improved plots of the drag on a sphere (Fig. [3]), by arguing that the leading 
divergence must be removed to better compare experiments and predictions (Fig. H|). This 
same criticism applies to plots of the drag on a cylinder. In the case of the cylinder, Co goes 
as R^ 1 (with logarithmic corrections) as R — > (Eqn. [T9l) . This means we ought to plot 
CdR/Ak. This function tends to zero as R — > 0, so it is not necessary to plot CdR/4tt — 1, 
as in the case of the sphere. Figure [7J shows both Finn's and Tritton's data re-plotted with 



7 Tritton does caution that his measurements may be negatively biased at higher Reynolds number (R > 30). 
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FIG. 6 (Color online) TVitton's measurements of the drag on a cylinder (jTrittonl . 



19591 1. 



the leading divergence removed. 



In 1965, K. O. L. F. Jayaweera (jJayaweera and Masonl . 119651 ) undertook drag mea- 
surements of the drag on very long (but finite) cylinders. At very low Reynolds number 
(R < 0.135), his data are available in tabular form. At higher Reynolds number, they had 
to be digitized. His data, plotted with the leading divergence removed, are also shown in 
Figure 

The agreement amongst these experiments is excellent. Henceforth, Finn's data will not 
be plotted, as it exhibits larger experimental variations, and is surpassed by the experiments 
of Jayaweera and Tritton. Jayaweera's data exhibit the least scatter, and may be slightly 
better than Tritton's. However, both experiments have comparable, large ratios of cylinder 
length to width (the principle source of experimental error), and there is no a priori reason 
to favor one experimental design over the other. We consider these two experiments to be 
equivalent for the purposes of evaluating theoretical predictions. 



As with th e sphere, there are numeric al calculations, including: Underwood (jUnderwoodl . 



19691 ). Son (ISon and Hanrattyl . Il969l ). Kawaguti (IKawaguti and Jainl . Il966l ). Dennis 
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FIG. 7 (Color online) Summary of me asurements of the drag on a cylinder (jFinnl . 



Javaweera and Masonl . 



1965; 



Tritton . 



1953 ). 



1953; 



Dennis and Shmshoni 



( lAllen and Southwell 



12fia), Thorn flThoml . Il933 l. Apelt flApeltl . Il96lh . and Allen 



19551 ) . Of these, most treat only a few Reynolds numbers, none 



of which are sufficiently s mall. Other s , such as Allen and Dennis, have had their results 

^j). The only applic able studies are Kawaguti 



196c 



subsequently questioned (lUnderwoodl . 
( IKawaguti and Jainl . Il966l ). and Underwood (lUnderwoodl . Il969l ). Kawaguti has a calcula- 
tion only for R = 0.5, and is omitted. Underwood's results are in principle important and 
useful, but are only available in a coarse plot, which cannot be digitized with sufficient accu- 
racy. Consequently, no numerical results will be used for evaluating analytical predictions. 

There are many different experimental and numerical drag coefficient measurements. We 
will subsequently use only the best as benchmarks for evaluating the performance of theo- 
retical predictions. In the case of the sphere, the experimental measure ments of Maxworthy 



(MaxworthyJ, Il965l ) as well as the numerica 



calc ulations of Dennis ( IDennis and Walkerl . 



19711 ) and Le Clair (iLe Clair and Hamieled . [1970J) all extend to sufficientl y small R and 
possess sufficient accuracy. For the cylinder th e experiments of both Tritton (ITrittonl . Il959l ) 



and Jayaweera (IJayaweera and Mason 



19651 ) are both excellent. Although they exhibit 



25 



small differences, we cannot judge either to be superior, and we will compare both with 
theoretical results. 



B. Theoretical history 

Since these problems were posed by Stokes in 1851, there have been many attempts to 
solve them. All of these methods involve approximations, which are not always rigorous (or 
even explicitly stated). There is also considerable historical confusion over contributions and 
attribution. 8 Here we review and summarize the substantial contributions to the literature, 
focusing on what approximations are used, in both deriving governing equations and in their 
subsequent solution. We discuss the validity and utility of important results. Finally, we 
emphasize methodological shortcomings and how they have been surmounted. 



1. Stokes and paradoxes 



In the first paper on the subject, Stokes approximated R = in Eqn. [6] and so 



resulting equation (a problem equivalent to solving Eqn. [12] with R = 0) ( IStokea . Il85ll ). 



ved the 



After applying the boundary conditions (Eqn. UM . his solution is given in terms of a stream 
function by Eqn. [281 

(28) 



2r z - 3r + 
4 V r 



By substituting ip(r, fi) into Eqns. HH EU and 



;) f 1 -" 2 ) 



famous result of Stokes, given by Eqn. 



C 



D 



6tt 
~R 



or by using Eqn. I2bj) . we reproduce the 



(29) 



Stokes also tackled the two dimensional cylinder problem in a similar fashion, but could 
not obtain a solution. The reason for his failure can be seen by setting R = in Eqn. [91 and 
attempting a direct solution. Enforcing the sin 9 angular dependence results in a solution of 
the form ip(r, 9) = (C*ir 3 + C^r lnr + C^r + C±/r) sin#. Here Cj are integration constants. 
No choice of Cj will meet the boundary conditions Eqn. ffTUl) . as this solution cannot match 



For an explanation of confusion over early work, see Lindgren (jLindgren 



19991). Proudman and Pearson 



( Proudman and Pearson L 19571) also b egin their article with an insightful, nuanced discussion, although 
there are some errors (|Lindgrenl . 1 1999T ) . 
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the uniform flow at large r. The best one can do is to set C\ = 0, resulting in a partial 
solution: 

ip(r,8) = C \2r\a.r-r + -J sin# (30) 

Nonetheless, this solution is not a description of fluid flow which is valid everywhere. More- 
over, due to the indeterminable constant C, Eqn. [HD] cannot be used to estimate the drag 
on the cylinder. 

A more elegant way to see that no solution m a y exis t is through dimensional analysis 



( jHappel and Brenner 



1973 



Landau and Lifschitz 



19991 ) . The force per unit length may 



only depend on the cylinder radius, fluid viscosity, fluid density, and uniform stream velocity. 
These quantities are given in Table IHIl with M denoting a unit of mass, T a unit of time, 
and L a unit of len gth. From these quantities, one may form two dimensionless groups 
(jBuckinghaml . 11914 ) : n = R = \uoo\a/u, U\ = F^ et /(pu\u 00 \). Buckingham's IT Theorem 



([Buckingham! . 



19141 ) then tells us that: 



F(R) 



(31) 



If we make the assumption that the problem does not depend on R, as Stokes did, then we 
obtain II i = const, whence 

F Net oc pz/juool (32) 

However, Eqn. [32] does not depend on the cylinder radius, a! This is physically absurd, 
and demonstrates that Stokes' assumptions cannot yield a solution. The explanation is 
that when we take the R — > limit in Eqn. [3TJ we made the incorrect assumption that 
F(R) tended toward a finite, non-zero limit. This is an e xample of incompl ete similarity, or 
similarity of the second kind (in the Reynolds number) ( BarenblattJ . 119961 ). Note that the 
problem of flow past a sphere involves force, not force per unit length, and therefore is not 
subject to the same analysis. 

Stokes incorrectly took this nonexistence of a solution to mean that steady-state flow past 
an infinite cylinder could not exist. This problem, which is know n as Stokes' paradox, has 
been shown to occur with any unbounded two-dimensional flow (iKrakowski and Charned . 
19531 ). But such flows really do exist, and this mathematical problem has since been resolved 
by the recognition of the existence of boundary layers. 

In 1888, Whitehead, attempted to find higher approximations for flow past a sphere, ones 
which would be valid for small but non-negligible Reynolds numbers ( [Whitehead! . 118881 ) . He 
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Quantity 


Description 


Dimensions 


-?*Net 


Net Force per Unit Length 


MT- 2 


V 


Kinematic Viscosity 


l?T v 


a 


Cylinder Radius 


L 


P 


Fluid Density 


ML' 3 


l^oo I 


The Uniform Stream Speed 


LT- 1 



TABLE III Dimensional analysis of Stokes' problem. 



used Stokes' solution (Eqn. [28]) to approximate viscous contributions (the LHS of Eqn. [T2l) . 
aiming to iteratively obtain higher approximations for the inertial terms. In principle, this 
approach can be repeated indefinitely, always using a linear governing equation to obtain 
higher order approximations. Unfortunately, Whitehead found that his next order solution 
could not meet all of the boundary condit ions (Eqn. [H]), because he could not match the 
uniform stream at infinity (IVan Dykd . ll975l ). These difficulties are analogous to the problems 
encountered in Stokes' analysis of the infinite cylinder. 

Whitehead's approach is equivalent to a perturbative expansion in the Reynolds 
number, an approach which is "never valid in problems of uniform streaming" 
(jProudman and Pearson! 119571 ). This mathematical difficulty is common to all three- 
dimensional uniform flow problems, and is known as Whitehead's paradox. Whitehead 
thought this was due to discontinuities in the flow field (a "dead-water wa ke"), but this 
is incorrect, and his "paradox" has also since been resolved (IVan Dykd . Il975l ). 



2. Oseen's equation 



a. Introduction 

In 1893, Rayleigh pointed out that Stokes' solution would be uniformly applicable if 
certain inertial forces were included, and noted that the ratio of those inertial forces to the 
viscous forces w hich Stokes considere d could be used to estimate the accuracy of Stokes' 
approximations (Lord Rayleigh . 1893 ). 

Building on these ideas in 1910, C. W. Oseen proposed an ad hoc approximation to the 
Navier-Stokes equations which resolved both paradoxes. His linearized equations (the Oseen 
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equations) attempted to deal with the fact that the equations governing Stokes' perturbative 
expansion are invalid at large \ f\, where they neglect important inertial terms. In addition 
to Oseen, a number of workers have applied his equations to a wide variety of problems, 
including both the cylinder and the sphere. 9 

Oseen's governing equation arises independently in several different contexts. Oseen 
derived the equation in an attempt to obtain an approximate equation which describes the 
flow everywhere. In modern terminology, he sought a governing equation whose solution is 
a uniformly valid approximation to the Navier-Stokes equations. Whether he succeeded is 
a matter of some debate. The short answer is "Yes, he succeeded, but he got lucky." 

This story is further complicated by historical confusion. Oseen's equations "are valid but 



for the wrong reason" (jLindgren . 



19991 ); Oseen originally objected to working in the inertial 



frame where the solid body is at rest, and therefore undertook calculations in the rest frame of 
uniform stream. This complication is overlooked largely because many subsequent workers 
have only understood Ose en's intricate three paper analysis through the lens of Lamb's 
later work (jLamb H.I . I1911I ). Lamb — in addition to writing in English — presents a clearer, 
"shorter way of arriving at his [ Oseen's] results" , which he characterizes as "somewhat long 



and intricate." (jLamb Hi Il91ll ) 



In 1913 Fritz Noethe r, using both R ayleigh's and Oseen's ideas, analyzed the problem 



using stream functions ((Noether 



19131 ). Noether's paper prompted criticisms from Oseen, 



who then revisited his own work. A few months late r, Oseen published another paper, which 
included a new result for Co (Eqn. [39]) f Osee J . 1913 ). Burgess also explains the development 
of Oseen's equation, and present s a clear derivation o f Oseen's principal results, particularly 
of Oseen's new formula for Cr> ( jBurgess R.W.I . Il916l ). 



Lindgren offers a detailed discussion of these historical developments (jLindgrenl . Il999f ). 
However, he incorrectly reports Noether's publication date as 1911, rather than 1913. As 
a result, he incorrectly concludes t hat Noether's work was independent of Oseen's, and 
contradicts claims made in Burgess (jBurgess R.W.I . Il916l ). 



Although the theoretical justification for Oseen's approximations is tenuous, its success at 



9 Lamb (jLambl . Il932f ) solved the Oseen equations for the cylinder approximately, as Oseen (Oseen, 1910) 
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7) 



did for the sphere. The Oseen eq uations have been solved exactly for a cylinder by Faxen (jFaxenl 
as well as by Tomotika and A oi (jTomotika and Aoil . |l95fj), and those for the sphere were solved exactly 
by Goldstein ( Goldsteinl . [l929 ). 
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resolving the paradoxes of both Stokes and Whitehead led to widespread use. Oseen's equa- 
tion has been fruitfully substituted for the Navier-Stokes' equations in a broad array of low 
Reynolds number problems. Happel and Brenner describe its applicat ion to many problems 
in th e dynamics of small particles where interactions can be neglected (IHappel and Brenner! . 



19731 ). Many workers have tried to explain the utility and unexpected accuracy of Oseen's 



governing equations. 

Finally, the Oseen equation, as a partial differential equation, arises in both matched 
asymptotic calculations and in our new work. In these cases, however, its genesis and 
interpretation is entirely different, and the similarity is purely formal. Due to its ubiquity 
and historical significance, we now discuss both Oseen's equation and its many different 
solutions in detail. 

b. Why Stokes' approximation breaks down 

Oseen solved the paradoxes of Stokes and Whitehead by using Rayleigh ' s insig ht: compare 
the magnitude of inertial and viscous forces (ILord Rayleighl . Il893l : lOseenl . Il910l ). Stokes and 
Whitehead had completely neglected inertial terms in the Navier-Stokes equations, working 
in the regime where the Reynolds number is insignificantly small (so-called "creeping flow"). 
However, this assumption can only be valid near the surface of the fixed body. It is never 
valid everywhere. 

To explain why, we follow here the s pirit of Lamb's analysis, presenting Oseen's conclu- 
sions "under a slightly different form." (jLamb H.I . Il91ll ) 

Consider first the case of the sphere. We can estimate the magnitude of the neglected 
inertial terms by using Stokes' solution (Eqn. UM . Substituting this result into the RHS of 
Eqn. [12], we see that the dominant inertial components are convective accelerations arising 
from the nonlinear terms in Eqn. [121 These terms reflect interactions between the uniform 
stream and the perturbations described by Eqn. [28] For large values of \r\, these terms are 
of 0{Rr- 2 ). 

Estimating the magnitude of the relevant viscous forces is somewhat trickier. If we 
substitute Eqn. [28] into the LHS of Eqn. [121 the LHS vanishes identically. To learn anything, 
we must consider the terms individually. There are two kinds of terms which arise far from 
the sphere. Firstly, there are components due solely to the uniform stream. These are 
of 0(r~ 2 ). However, the uniform stream satisfies Eqn. [T2l independently, without the 
new contributions in Stokes' solution. Mathematically, this means that all of the terms of 
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O (r~ 2 ) necessarily cancel amongst themselves. 10 We are interested in the magnitude of the 
remaining terms, perturbations which result from the other components of Stokes' solution. 
These viscous terms (i.e. the d^j term in Eqn. [T^i) are of O (r -3 ) as r — > oo. 

Combining these two results, the ratio of inertial to viscous terms, in the r — > oo limit, 
is given by Eqn. I3"3"l 

inertial _ . „ . . , 

= 0(Rr) 33 

viscous 

This ratio is small near the body (r is small) and justifies neglecting inertial terms in 
that regime. However, Stokes' implicit assumption that inertial terms are everywhere small 
compared to viscous terms breaks down when Rr ~ O ( 1 ) , and the two kinds of forces are 
of the same magnitude. In this regime, Stokes' solution is not valid, and therefore cannot be 
used to estimate the inertial terms (as Whitehead had done). Technically speaking, Stokes' 
approximations breaks down because of a singularity at infinity, an indication that this is a 
singular perturbation in the Reynolds' number. As Oseen pointed out, this is the genesis of 
Whitehead's "paradox" . 

What does this analysis tell us about the utility of Stokes' solution? Different opinions 
can be found in the literatu re. Happel, for instance, claims that it "is not uniformly valid" 
(IHappel and Brennerl . 119731 1 . while Proudman asserts "Stoke s' solution is therefore actually 



a uniform approximation to the total velocity distribution." (IProudman and Pearson! . 119571 ) 
By a uniform approximation, we mean that the approxim ation asymptotically app r oache s 
the exact solution as the Reynolds' number goes to zero (IKaplun and Lagerstroml . 119571 ); 
see Section III. CI for further discussion. 

Proudman and Pearson clarify their comment by noting that although Stokes' solution is a 
uniform approximation to the total velocity distribution, it does not adequately characterize 
the perturbation to the uniform stream, or the derivatives of the velocity. This is a salient 
point, for the calculations leading to Eqn. [33] examine components of the Navier-Stokes 
equations, not the velocity field itself. These components are forces — derivatives of velocity. 

However, Proudman and Pearson offer no proof that Stokes' solution is actually a uniform 
approximation, and their claim t hat it is "a valid approximation to many bulk properties of 



the flow, such as the resistance" (IProudman and Pearson 



19571 ) goes unsupported. In fact 



10 



VanDvke (Van Dvke 




1975) 


( Proudman and Pearsor 


l. 1957) 
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any calculation of the drag necessitates utilizing derivatives of the velocity field, so their 
argument is inconsistent. 

We are forced to conclude that Stokes' solution is not a uniformly valid approximation, 
and that his celebrated result, Eqn. [2H] is the fortuitous result of uncontrolled approxima- 
tions. Remarkably, Stokes' drag formula is in fact the correct zeroth order approximation, 
as can be shown using either matched asymptotics or the Oseen equation! This coincidence 
is essentially due to the fact that the drag is determined by the velocity field and its deriva- 
tives at the surface of the sphere, where r = 1, and Eqn. [33] is O (R 1 )- The drag coefficient 
calculation uses Stokes' solution in the regime where his assumptions are the most valid. 

A similar analysis affords insight into the origin of Stokes' paradox in the problem of 
the cylinder. Although we have seen previously that Stokes' approach must fail for both 
algebraic and dimensional considerations, examining the ratio between inertial and viscous 
forces highlights the physical inconsistencies in his assumptions. 

We can use the incomplete solution given by Eqn. [30]to estimate the relative contributions 
of inertial and viscous forces in Eqn. [9j More specifically, we examine the behavior of these 
forces at large values of r. Substituting Eqn. [30] into the RHS of Eqn. EJ we find that the 
inertial forces are O (RC 2 log r/r 2 ) as r — > oo. 

We estimate the viscous forces as in the case of the sphere, again ignoring contributions 
due solely to the uniform stream. The result is that the viscous forces are O (Clog r/r 3 ). 11 
Combining the two estimates, we obtain the result given in Eqn. [341 

inertial . „ . , . 

0(Rr) (34) 



VISCOUS 

This result demonstrates that the paradoxes of Stokes and Whitehead are the result of 
the same failures in Stokes' uncontrolled approximation. Far from the solid body, there is a 
regime where it is incorrect to assume that the inertial terms are negligible in comparison 
to viscous terms. Although these approximations happened to lead to a solution in the case 
of the sphere, Stokes' approach is invalid and technically inconsistent in both problems. 

c. How Oseen Resolved the Paradoxes 



ii 



This result disag rees with the results of Proudman (jProudman and Pearsonl I1957T ) and VanDyke 



( Van Dvkel . Il97a ). who calculate that the ratio of inertial to viscous forces ~ Rr In r. However, both 



results lead to the same conclusions. 
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Not only did Oseen identify the physic al origin for t he breakdowns in previous approxi 



mations, but he also discovered a solution (lOseen 



19101 ) . As explained above, the problems 



arise far from the solid body, when inertial terms are no longer negligible. However, in this 
region (r ^> 1), the flow field is nearly a uniform stream — it is almost unperturbed by the 
solid body. Oseen's inspiration was to replace the inertial terms with linearized approxi- 
mations far from the body. Mathematically, the fluid velocity u in Eqn. [6] is replaced by 
the quantity Uoo + u, where u represents the perturbation to the uniform stream, and is 
considered to be small. Neglecting terms of O (|w| 2 ), the viscous forces of the Navier-Stokes' 
equation — R {u • Vu) — are approximated by R {uoo ■ Vu). 
This results in Oseen's equation: 

R (woo • W) = - Vp + V 2 u (35) 

The lefthand side of this equation is negligible in the region where Stokes' solution applies. 
One way to see this is by explicitly substituting Eqn. [28] or Eqn. [30] into the LHS of Eqn. 
[351 The result is of O (R). This can also be done self-consistently with any of the solutions 
of Eqn. [35] it can thereby be explicitly shown that the LHS can only becomes important 
when r> 1, and the ratios in Eqns. [331 and [3^1 are of O (1). 

Coupled with the continuity equation (Eqn. [5]), and the usual boundary conditions, the 
Oseen equation determines the flow field everywhere. The beautiful thing about Oseen's 
equation is that it is linear, and consequently is solvable in a wide range of geometries. In 
terms of stream functions, the Oseen equation for a sphere takes on the form given by Eqn. 
The boundary conditions for this equation are still given by Eqn. [1] 



D^ = R[ + ) D^{r^) (36) 



r d/j dr / 
Here, D is defined as in Eqn. [12] 

For the cylinder, where the boundary conditions are given by Eqn. [TO], Oseen's equation 
takes the form given by Eqn. 



Vty(r, 9) = R (cos(0)|; - V^(r, 9) (37) 

Here V is defined as in Eqn. [9] This equation takes on a particularly simple form in 
Cartesian coordinates (where x = rcos6): (V 2 — Rd x ) V 2 ip(r, 9) = 0. 

A few historical remarks must be made. First, Oseen and Noether were motivated to 
refine Stokes' work and include inertial terms because they objected to the analysis being 
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done in the rest frame of the solid body. While their conclusions are valid, there is nothing 
wrong with solving the problem in any inertial frame. Secondly, Oseen made no use of 
stream functions; the above equations summarize results from several workers, particularly 
Lamb. 

There are many solutions to Oseen's equations, applying to different geometries and 
configurations, including some exact solutions. However, for any useful calculations, such 
as Cfl, even the exact solutions need to be compromised with approximations. There have 
been many years of discussion about how to properly interpret Oseen's approximations, and 
how to understand the limitations of both his approach and concomitant solutions. Before 
embarking on this analysis, we summarize the important solutions to Eqns. EH and 

d. A plethora of solutions 

Osee n himself pr ovided the first solution to Eqn. [36J, solving it exactly for flow past a 
sphere ( lOseenl . 1 1 9 1 Of ) . Eqn. [381 rep roduces this result in terms of stream functions, a formula 
first given by Lamb (jLambl . Il932l ). 



Mr,6) = ~ ( 2r 2 + - ] sin 2 #- — (l + cos#) (l 
4 \ r J 2R \ 



e 2 



ii?r(l-cos0) 



(38) 



This solution is reasonably behaved everywhere, and may be used to obtain Oseen's improved 
approximation for the drag coefficient (Eqn. I3"9|) . 

3 



C 



D 



6tt 



1 + 1r\ + O (R 2 ) 



(39) 



Oseen obtained t his prediction for Cd after the prompting of Noeth er, and only pr e sente d 



it in a later paper (jOseenj, Il913l ). Burgess also obtained this result (IBurgess R.W.I . Il916l ). 
Oseen's work was hailed as a resolution to Whitehead's paradox. While it did resolve the 
paradoxes (e.g., he explained how to deal with inertial terms), and his solution is uniformly 
valid, it does not posses sufficient accuracy to justify the "3/8/2" term in Eqn. [391 What 
Oseen really did was to rigorously derive the leading order term, proving the validity of 
Stokes' result (Eqn. I29"j) . Remarkably, his new term is also correct! This is a coincidence 
which will be carefully considered later. 

This solution (Eqn. [38]) is exact in the sense that it satisfies Eqn. [361 However, it does 
not exactly meet the boundary conditions (Eqn. [T4"|) at the surface of the sphere. It satisfies 
those requirements only approximately, to 0(R l ). This can readily be seen by expanding 
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Eqn. [38] about r = 1: 



J ^2r 2 - 3r + sin 2 + O (R 1 ) 



(40) 



Up to O (R) this is simply Stokes' solution (Eqn. 128]) . which vanishes identically at r = 1. 
The new terms fail to satisfy the boundary conditions at the surface, but are higher order in 
R. Thus Oseen's solution is an exact solution to an approximate governing equation which 
satisfies boundary conditions approximately. The implications of this confounding hierarchy 
of approximations will be discussed below. 

Lamb contrib uted a simplified method for both deriving and solving Ose en's equation 



(lLamb H. 



Goldstein 



1911 ). His formulation was fru itfully used by later workers (e.g., (jFaxenl . 1 192 7 ; 



1929; 



Tomotika and Aoi 



19501 ) ). and Lamb himself used it both to reproduce 



Oseen's results and to obtain the first result for the drag on an infinite cylinder. 

Lamb's basic solution for flow around an infinite cylinder appears in a number of guises. 
His original solution was given in terms of velocity components, and relied on expansions 
of modified Bessel functions which kept only the most important terms in the series. This 
truncation results in a solution (Eqn. |4T]) which only approximately satisfies the governing 
equations (Eqn. [57]) . and is only valid near the surface. 



Uy 



1 +5 ( 7 - - + log— + - (r -1]— logr 



dx 2 



6 
2 




(r 2 - 1) 



d 2 



dxdy 



logr 



(41a) 

(41b) 
(41c) 



In this equation, 5 = f~ — 7 — log j) . 

Note that, although it only approximately satisfies Oseen's governing equation, this result 
satisfies the boundary conditions (Eqn. H]) exactly. Lamb used his solution to derive the 
first result (Eqn. 142]) for the drag on an infinite cylinder, ending Stokes' paradox: 

47T 



c 



D 



R 



(6) 



(42) 



In his own words, "... Stokes was led to the conclusion that steady motion is impossible. It 
will appear that when the inertia term s are part i ally ta ken into account ... that a definite 



value for the resistance is obtained." (lLamb H. 



19111 ) As with all analysis based on the 



ad-hoc Oseen equation, it is difficult to quantify either the accuracy or the limitations of 
Lamb's result. 



35 



Many authors formulate alternate expressions of Lamb's solution by retaining the mod- 
ified Bessel functions rather than replacing them with expansions valid for small R and r. 
This form is giv e n by Eqn. |4"3"| and is related to the incomplete form given by VanDyke (p. 



162) flVan Dvkel . Il975f ) 



12 



U , 



u, 



u, 



, x 1 2x 
i + S[^-—; + —, 



r / xy 




2y_ 

Rr 2 



-7 efe/2 ^(f)) (43a) 

(£)) (43b) 



(43c) 



Here I n and K n are modified Bessel functions. 

In contrast to Eqn. HU this solution is an exact solution to Oseen's equation (Eqn. [3 
but only meets the boundary conditions to first approximation. In particular, it breaks 
down for harmonics other than sin 9. Whether Eqn. [41] or Eqn. 03] is preferred is a matter 
of some debate, and ultimately depends on the problem one is trying to solve. 

Some workers prefer expressions like Eqn. H3], which are written in terms of u. Unlike the 
solutions for the stream function, these results can be written in closed form. This motivation 
is somewhat misguided, as applying the boundary conditions nonetheless requires a series 
expansion. 

In terms of stream functions Eqn. H3] transforms into Eqn. HH (IProudman and Pearson! . 



19571 ). 



5 

r H | sin I 

2r 



n=l 



Rr\ rsmnO 
~Y) n 



(44) 



Here, 



4> n (x) = 2K x {x)I n (x) + K (x) (I n +i{x) + I n -i{x)) 
This result is most easily derived as a special case of Tomotika's general solution (Eqn. 



9jj (ITomotika and Aoi 



1950) , although Proudman et al. intimate that it can also be directly 



derived from Lamb's solution (Eqn. |4*3|) (IProudman and Pearson! . 1 195 71 ). 

Bairstow et al. w ere the first to retain Bessel fu nctions while solving Oseen's Eqn. for 



flow past a cylinder (jBairstow and Cave B.M.I . Il923l ). They followed Lamb's approach, but 



endeavored to extend it to larger Reynolds' numbers, and obtained the drag coefficient given 



Note that VanDyke incorrectly attributes to this result to Oseen, rather than to Lamb. 
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in Eqn. H5] When expanded near R = 0, this solution reproduces Lamb's result for Co 
(Eqn. H21 . It can also be obtained from Tomotika's more general solution (Eqn. H9|) . 

° D = R(I (R/2)K (R/2) + htRfflKjRffl) (45) 
Bairstow also made ext ensive com parisons between experimental measurements of Co and 



theoretical predictions (jReli Il914l ). He concluded, "For the moment it would appear that 
the maximum use has been made of Oseen's approximation to the equations of viscous fluid 
motion." 

At this point, the "paradoxes" were "resolved" but by an approximate governing equation 
which had been solved approximately. This unsatisfactory state of affairs was summarized 
by Lamb in the last edition of his book: "... even if we accept th e equa t ions a s adequate the 
boundary-conditions have only been approximately satisfied." (iLambl . 119321 ) His comment 
was prompted largely by the work of Hilding Faxen, who initiated the next theoretical 
development, exact solutions to Oseen's approximate governing equation (Eqn. [35]) which 
also exactly satisfy the boundary conditions. 

Beginning with his thesis and spanning a number of papers Faxen systemat ically investi 



gated the application of boundary conditions to solutions of Oseen's equations (jFaxenl . 



1921 



19231 ). Faxen initially studied low Reynolds number flow around a sphere, and he began 



by re-examining Oseen's analysis. He derived a formula for Co which differed from Oseen's 
accepted result (Eqn. [39]). Faxen realized that this was due to differences in the application 
of approximate boundary conditions; within the limitations of their respective analyses, the 
results actually agreed. 

Faxen next solved Oseen's equation (Eqn. [36]) . but in bounded, finite spaces where the 
boundary conditions could be satisfied exactly. He initially studied flow near infinite parallel 
planes, but ultimately focused on flow around a sphere within a cylinder of finite radius. 
He aimed to calculate the drag force in a finite geometry, and then take the limit of that 
solution as the radius of the cylinder tends to infinity. 

Unfortunately, in the low Reynolds number limit, the problem involves incomplete simi- 
larity, and it is incorrect to assume that solutions will be well behaved (e.g., tend to a finite 
value) as the boundary conditions are moved to infinity. 

The drag force which Faxen calculated involved a number of undetermined coefficients, 
so he also calculated it using solutions to Stokes' governing equations. This solution also has 
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unknown coefficients, which he then calculated numerically. Arguing that the two solutions 
ought to be the same, he matched coefficients between the two results, substituted the 
numerical coefficients, and thereby arrived at a drag force based on the Oseen governing 
equation. 

This work is noteworthy for two reasons. First, the matching of coefficients between 
solutions derived from the two different governing equations is prescient, foreshadowing the 
development of matched asymptotics 30 years later. Secondly, Faxen ultimately concluded 
that O s een's "improvement" (Eqn. [39]) on Stokes' drag coefficient (Eqn. [29!) is invalid 
(jFaxenl . Il923l ). Faxen's analysis demonstrates that — when properly solved — Oseen's 



equat ion yields the same drag coefficient as Stokes', without any additional terms (ILindgrenl . 



19991). 



Studies by Bohlin and 



Haberman and Saure . 



1958 



iaberman conc ur with Faxen's conclusions (jBohlinl . 



1960 ; 



Lindgrenl . Il999h . It is not surprising that his results reject 



Oseen's new term (3R/8). We previously explained that Oseen's analysis, although it elim- 
inates the "paradoxes", does not posses sufficient accuracy to justify more than the lowest 
order term in Eqn. 1391 

However, Faxen's results suffer from problems. First, they cannot be systematically used 
to obtain better approximations. Secondly, Faxen actually solves the problem for bounded 
flow, with the boundary conditions prescribed by finite geometries. He uses a limiting 
procedure to extend his solutions to unbounded flow (with boundary conditions imposed 
on the uniform stream only at infinity, as in Eqn. H|). In problems like this, which involve 
incomplete similarity, it is preferable to work directly in the infinite domain. 

Faxen's meticulous devotion to properly applying boundary conditions culminated in 
the first complete solution to Eqn. [371 In 1927, he published a general solution f or flow 
aroun d an infinite cylinder which could exactly satisfy arbitrary boundary conditions (jFaxenl . 



19271 ). Unfortunately, Faxen's solution contains an infinite number of undetermined integra- 



tion constants, and approximations must be used to determine these constants. Although 
this destroys the "exact" nature of the solution, these approximations can be made in a 
controlled, systematic fashion — an improvement over the earlier results of Lamb and Os- 
een. Although Faxen's heroic solution was the first of its kind, his real insight was realizing 
that approximations in the application of boundary conditions could be as important as the 
approximations in the governing equations. 
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His formal solutions are in essence a difficult extension of Lamb's reformulation of Os- 
een's equations, and they inspired several similar solutions. In 1929, Goldstein completed 



a similarly hercu 



around a sphere (jGoldsteinl . 



ean calcu l ation to derive a general solution to Oseen's equation for flow 



1929 ). Like Faxen's result for the cylinder, Goldstein's solution 



can — in principle — exactly satisfy the boundary conditions. Unfortunately, it also suffers 
from the same problems: It is impossible to determine all of the infinitely many integration 
constants. 

Goldstein's soluti on is summarized by Tom otika, who also translated it into the language 
of stream functions (ITomotika and Aoil . Il950l ). We combine elements from both papers in 
quoting the solution given in Eqn. 



V(r,0) = -r 2 Qa(cos0) + ^ ( B n r- n + ^ X m r 2 $ m , n (rR/2) ) Q n (cosl 
In this equation, 



n=l 



m=0 



Qn(fJ-) 



(46) 



(47a) 
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Xm-l{X) + ~ —Xm+l{x) fm,n{ X ) 



2m - 1 

m 



2m + 3 
m + 1 



2m + l /m - 1,n(x) + 2/;; ] 



f m+1) n(x) Xm{x) 



fm,n (*£} 



(2m + 1) 



\2xJ 



(47b) 
(47c) 



(2n + l)J2 



(2j)\{2m-2j)\(2n-2j)\ 



x 



j=0 

(m + n — j 



(m-j)\(n-j)\ 



(j!) 2 (2m + 2n - 2j + 1)! 

2 

4>m+n-2j(x) 



(47d) 



<Pn{x) 
fm,n (•£-) 



^ ] C m (k) 



(47e) 



Here i^ n (^) and J n (a;) are Bessel functions, P m (x) are Legendre polynomials, and C m (k) is 
the coefficient of x k in P m (x). Note that the s econd expression for / mj „(x), written in terms 
of derivatives, is computationally convenient (jGoldsteinl . Il929i ). 

Eqn. HB] is given with undetermined constants of integrat ion, B n and X m . Metho ds to 
determine these cons t ants w ere discussed by both Tomotika (ITomotika and Aoil . Il950l ) and 



Goldstein (IGoldstein 



19291 ) . We will present our own analysis later. 
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There are many different results which have been obtained using the above general so- 
lution. The exact formula for the stream function and the drag coefficient depend on what 
terms in the solution are retained, and how one meets the boundary conditions. In general, 
retaining n angular terms in Eqn. HH1 requires the retention of m = n — 1 terms in the 
second sum. In his original paper, Goldstein retains three terms in each series, and thereby 
calculates the formula for Cd given in Eqn. HHJ 
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8 320 2560 



»i™ :fl . +i ^w fl . +0(Jf)1 (48) 



2150400 



17203200 



The coefficient of the last term reflects a correction due to Shanks (Shanks 



19551 ). 



To obtain the result in Eqn. HHJ Goldstein both truncated his solution for the stream 
function and then expanded the resulting Cd about R = 0. Van Dyke extended this result 
to include an additional 24 terms, for purposes of studying t he mathematical structure of 



the series, but not because of any intrinsic physical meaning ([Van Dykd . Il970l ). Van Dyke 
does not state whether he was including more harmonics in the stream function solution or 
simply increasing the length of the power series given in Eqn. 



In addition to expressing Goldstein's solution for the geometry of a sphere in terms of 



stream functions, T o moti 



and Goldstein ( 



sa derived his own exact solution to Eqn. [37] for flo w past a cy 



( ITomotika and Aoil. 11950 ) . Tom otika closely followed the spirit of Lamb (ILamb H 



Goldstein 



inder 



19 111 ) 



19291 ). and his resulting "analysis is quite different from Faxen's." 



( ITomotika and Aoil . Il950l ). His solution to Eqn. [37] is given in Eqn. [39] below, conveniently 
expressed in terms of stream functions. Note that Tomotika's result suffers from the same 
problems as his predecessors: An infinite number of undetermined integration constants. 



ip(r,6) = rsinO + I B„ r - n + X m r$ m , n (rR/2) I sinnO 



(49) 



n=l 



m=0 



Where 



$m,n( X ) = (K m+1 (x) + K m ^x{x)) (J TO _ n (z) + I m+n (x)) 

As before, B n and X m are constants of integration which need to be determined by the 
boundary conditions (Eqn. [TUT) . 

As with Goldstein's solution for the sphere, approximations are necessary in order to 
actually calculate a drag coefficient. By retaining the m = and n — 1 terms, Tomotika 



40 



reproduced Bairstow's result for Cd (Eqn. 1451) . He also numerically calculated drag coef- 
ficients based on retaining more terms. As with the Goldstein solution, keeping n angular 
terms requires keeping m = n — 1 terms in the second sum. 

The solutions given in Eqns. US] and HH] represent the culmination of years of efforts to 
solve Oseen's equation for both the sphere and the cylinder. These general solutions are 
also needed in both matched a symptotics and the new techniques presented in this section 



(jProudman and Pearson! . 



19571 ) 



There is a final noteworthy solution to Eqn. [37J In 1954, Imai pu blished a genera l 
method for solving the problem of flow past an arbitrary cylindrical body (jlmai Isad . 1 19541 ). 
His elegant technique, based on analytic functions, applies to more general geometries. Imai 
calculated a formula for Cd, approximating the functions in his exact solution with power 
series about R = 0. His result (re-expressed in our notation) is given in Eqn. [5TJ 

C D = — 5 + R[ 1 (51) 

R V 2 4 32 J K ' 

Note that Imai's result agrees with Eqn. I4"2l at lowest order, the only order to which Oseen's 

equation really applies. A priori, his result is neither better nor worse than any other solution 

of Oseen's equation. It is simply different. 

e. Discussion 

We have presented Oseen's governing equations for low Reynold number fluid flow. These 
equations are a linearized approximation to the Navier- Stokes' equations. We have also 
presented a number of different solutions, for both stream functions and drag coefficients; 
each of these solutions comes from a unique set of approximations. The approximations 
which have been made can be put into the following broad categories: 

• The governing equation — Oseen's equation approximates the Navier-Stokes equa- 
tions. 

• Solutions which only satisfy the Oseen's equation approximately. 

• Solutions which only satisfy the boundary conditions approximately. 

• Solutions where the stream function is expanded in a power series about R = after 
its derivation. 

• Approximations in the drag coefficient derivation. 
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• Drag coefficients which were expanded in a power series about R = after their 
derivation. 

The first approximation is in the governing equations. Oseen's approximation is an ad 
hoc approximation which, although it can be shown to be self-consistent, requires unusual 
cleverness to obtain. Because it is not derived systematically, it can be difficult to understand 
either its applicability or the limitations of its solutions. There have been years of discussion 
and confusion about both the equation and its solutions. The short answer is this: Oseen's 
governing equation is a zeroth order uniformly valid approximation to the Navier Stokes 
equation; the equation and its solutions are valid only a t O (R °) . 



It is not easy to prove this claim rigorously (iFaxenl . Il923l ). However, it can be easily 
shown that Oseen's equations are self-consistent with its solutions, and that the error in 
the solution is of 0(R l ). One way to explicitly demonstrate this is by substituting a 
solution of Oseen's equation into the LHS of the Navier-Stokes equations (Eqn. [6]), thereby 
estimating the contribution of inertial terms for the flow field characterized by the solution. 
By repeating that substitution into the LHS of Oseen's equation (Eqn. [35]), one can estimate 
the contribution of inertial terms under Oseen's approximations. Comparing the two results 
gives an estimate of the inaccuracies in Oseen's governing equations. 

Concretely, for the sphere, we substitute Eqn. [38] into the RHS of Eqn. [361 an d into the 
RHS of Eqn. [12] The difference between the two results is of O {R 1 )- 

For the cylinder, substitute Eqn. [44] into the RHS of Eqns. [37] and [9] The difference 
between the exact and approximate inertial terms is of O (R5), where 5 is defined as in Eqn. 
US 

These conclusions do not depend on the choice of solution (or on the number of terms 
retained in Eqn. |4"4"|) . They explicitly show that the governing equation is only valid to O (R) 
(or O (RS)). Consequently, the solutions can only be meaningful to the same order, and the 
boundary conditions need only be satisfied to that order. With these considerations, almost 
all of the solutions in the preceding section are equivalent. The only ones which are not - 
such as Eqn. [41] — are those in which the solution itself has been further approximated. 13 

Since the formulae for determining Co (Eqns. [T9l and [241) are of the form 1/R + terms 



13 In this case, the Bessel functions have been expanded near R = and are no longer well behaved as 
R — > oo. 
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linear in stream function + nonlinear terms, a stream function which is valid to O (R) will 
result in a drag coefficient which is valid to O (R°). Thus, in all of the formula for Cd which 
have been presented so far, only the first term is meaningful. For a sphere, this is the Stokes' 
drag (Eqn. |2"§]K and for the cylinder, Lamb's results (Eqn. |4"2"]) . 

We have concluded that it is only good fortune that Oseen's new "3/8i?" term is actually 
correct. This concurs with the analysis of Proudman et al., who wrote, "Strictly, Oseen's 
method gives only the leading term ... and is sca rcely to be counted as superior to Stokes's 
method for the purpose of obtaining the drag." (IProudman and Pearson! . 119571 ) Proudman 
and Pearson also note that the vast effort expended finding exact solutions to Oseen's equa- 
tion is "of limited value." Goldstein's formula for Cd, for instance, is expanded to O {R 5 ), 
well beyond the accuracy of the original governing equations. The reason for Oseen's good 
fortune is rooted in the symmetry of the problem. Chester and Van Dyke both observe that 
the non-linear terms which Oseen's calculation neglects, while n eeded fo r a correct stream 



funct ion, do not contribute to the drag because of symmetry (jChesterl . 



1962 



Van Dykd . 



1975|). 



Lindgren argues that Faxen proved that, when the boundary conditions are met properly 
and Oseen ' s equ ations solved exactly, the resulting Cp is that obtained by Stokes (Eqn. 129]) 
(ILindgrenl . Il999l ). Whether this argument is correct does not matter, as Oseen's additional 
term is beyond the accuracy of his governing equations. 

There is another approximat ion which arises wh i le com puting; Cd in the context of Oseen's 
equation. Many workers (e.g., (ITomotika and Aoil . Il950l )) compute the pressure in Eqns. [HI 
and [21] by integrating Oseen's equation (Eqn. [35] rather than the Navier-Stokes equations 
(Eqn. [6]). In Eqns. [20] and [25] we presented a pressure calculation based on the Navier 
Stokes equations. Calculating pressure using the linearized Oseen equation introduces an 
additional approximation into Co- While not necessarily problematic or inconsistent, this 
approximation can be difficult to identify. 

/. Two different interpretations One criticism of the Oseen equation is that it may be 
obtained by linearizing the Navier-Stokes equations, without regard to the magnitude of 
inertial and viscous terms. By writing u = + 5u, treating 5u as a small perturbation, and 
expanding Eqn. [6] one can formally reproduce Oseen's equations. Clearly, the disturbance 
to the uniform stream is not negligible near the surface of the solid body, and therefore 
Oseen's equations "would appear to be a poor approximation in the neighborhood of the 
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body where the boundary condition u = requires that the true inertial term be small." 



(IHappel and Brenner 



1973|) 



This incorrect argument, put forth as a reason to use Stokes' solutions, overlooks the 
origins of Oseen's equations. The point of Oseen's approximation is that inertial terms 
are only significant at large values of \r\, where R\r\ is no longer negligible. Near the 
surface of the solid, the approximate inertial terms which Oseen introduced are negligible in 
comparison to the viscous terms, because they are multiplied by the factor R (in the LHS of 
Eqn. 1351) . Hence the difference between Oseen's and Stokes' equations in the neighborhood 
of the sphere will be of O (R), and is beyond the scope of either theory. 

g. Better approximations 

The approach of W hitehead was esse ntially to improve Stokes' solution for the sphere 
in an iterative fashion (jWhiteheadl . Il888l ). By substituting the first approximation into the 
governing equations, he estimated the neglected terms. He then tried, and failed, to solve 
the resulting governing equation. This approach fails because the Stokes' equations are not 
uniformly valid to zeroth order. 

Oseen's equations are uniformly valid, and, as Proudman remarked, "there seems little 
reason to doubt that Whitehead's iterative method, using Oseen's equation rather than 
Stokes's equation would yield an expansion, each successive term of which would represent a 
uniformly valid higher approximation to the flow. In each step of the iteration, a lower-order 
approximation would be used to calculate those particular inertia terms that are neglected 
... the exp ansion generated in thi s way would seem to be the most economic expansion 
possible." (jProudman and Pearson! . 119571 ) 

Proudman did not follow through on this idea, instead developing a solution based on 
matched asymptotics expansions (see below). In an app endix, Van Dyke relates the unpub- 
lished work of C. R. Illingworth (1947) ( VanDyke, 1975^ . Illingworth carried through White- 
head's program, deriving a new expression (Eqn. l52l) for Cd, which agrees to 0(R 2 lnR) 
with the later results of matched asymptotic calculations (Eqn. 1551) . 



C 



D 



6tt 
~R 



(} + l R+ i) R2 hgR + °- 1333i?2 + ^ R31o & R ~ °- 0034R3 + ■■■) ( 52 ) 



Although this result has since been subsumed by matched asymptotics, it is nonetheless 
remarkable, substantially improving any previous drag calculations, and rigorously justifying 
Oseen's 3/8R term. 
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There have also been efforts (e.g., (jShanksl . Il955l ; IVan Dykd . Il970l )) to "re-sum" Gold- 
stein's series expansion for Cd (Eqn. US]) . However, these results have little intrinsic (as 
opposed to methodological) value, as Goldstein's result is only valid to O (R). If applied to 
more accurate approximations, such as Eqn. these methods could be worthwhile. Alas, 
even improved approximations lack a sufficient numbers of terms in the expression for Cp 
to make this practicable. 

h. Summary 

Simply put, Oseen's equations resolved the paradoxes of Stokes and Whitehead, put 
Stokes' results on firm theoretical ground, and led to the first solution for the drag on a 
cylinder. Although the Oseen equations happen to provide a uniformly valid first approxi- 
mation, it is difficult to extend this work to higher order approximations. 

Figure [8] compares the "predictions" of Oseen theory to experimental and numerical 
data for the drag on a sphere. Again, Oseen's first order theory is, strictly speaking, not 
adequate to make the predictions with which it is traditionally credited. The theoretical 
drag coefficients are roughly valid for R < 0.2, with Goldstein's solution (Eqn. 1481) being 
slightly better than Oseen's prediction (Eqn. [39]). All are clearly superior to Stokes' formula 
(Eqn. [29]). 

Figure [8] also shows the prediction of Illingworth's second order Oseen theory (Eqn. 1521) . 
Not surprisingly, it gives the best prediction of Co, particularly when compared to Dennis' 
numerical results. 

Figure [9] shows the important predictions of Oseen theory for the drag on an infinite 
cylinder. As with the sphere, the theory is only truly entitled to predict the lowest order 
term. Figure [9] shows decent agreement with the data. Although more "exact" solutions 
(such as Bairstow's and Imai's) do better than Lamb's lowest order solution, this is purely 
coincidental. Tomotika's solutions exhibit similar characteristics to these two solutions. 



3. Matched asymptotics 

Efforts to systematically improve Oseen's results led to the development of matched 
asymptotic expansions. 14 This branch of applied mathematics was developed gradually, 

14 This technique is also known as the method of inner and outer expansions or double asymptotic expansions. 
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Le Clair and Hamieled . 



Drag 



o n a sphere, e xperi ment vs. Oseen theory (jDennis and Walker 



1970; 



Maxworthvl . 



1971 



19651 ). The Stokes' solution (Eqn. [29)) is shown at the 



bottom for reference. In these coordinates, it is defined by the line y = 0. 



with systematic work beginning with papers by Kaplun and Lagerstrom et al. (IKaplunl . 



1954 



Lagerstrom and Cold . Il955l ). Kaplun subsequently used the se techniques to ca 



the drag on a cylinder, obtaining an entirely new result for Co ( ILagerstrom et al. 



culate 



1967L 



Proudman and Pearson subsequently applied matched asy mptotics to both the sphere an d 



the cylinder, deriving a new result for the drag on a sphere (jProudman and Pearson! . 



19571 ) 



"The principle of asymptotic matching is simple. The interval on which a 
boundary- value problem is posed is broken into a sequence of two or more over- 
lapping subintervals. Then, on each subinterval perturbation theory is used to 
obtain an asymptotic approximation to the solution of the differential equation 
valid on that interval. Finally, the matching is done by requiring that the asymp- 
totic approximations have the same functional form on the overlap of every pair 
or intervals. This gives a sequence of asymptotic approximations ... the end 
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result is an app roximate solution to a b oundary-value problem valid over the 



entire interval" (IBender and Orzag 



1999|). 



Both of the two low Reynolds number problems are attacked in similar fashion. The 
problem is divided into only two regions. The first region is near the surface of the solid 
body. In this region, inertial terms are small, the approximation of Stokes (R ~ 0) applies, 
and the problem is solved perturbatively (in R). At each order in R, the two no-slip boundary 
conditions at the surface are applied. One undetermined constant remains (at each order 
in R). Loosely speaking, it is determined by the boundary condition as \r\ — > oo. This 
expansion is referred to as the Stokes expansion. 

The second region is far from the sphere, where inertial terms are important. In this 
region, R\r\ ~ 0(1), and the approximations which led to Oseen's governing equation 
apply. The Oseen problem is then solved perturbatively, and the boundary condition as 
\r\ —>■ oo is applied. There are two undetermined constants remaining; they are related to 
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the boundary conditions on the surface. This perturbative expansion is referred to as the 
Oseen expansion. 

The next part of this calculation is asymptotic matching, which determines the remaining 
coefficients. 15 In this process, we expand the Oseen expansion for small R\r\, and the 
Stokes expansion for large \r\. By choosing the three hitherto undetermined coefficients 
appropriately, these two limiting forms are made to agree order by order in R. For this 
to be possible, the two asymptotic functional forms must overlap. With the coefficients 
determined, the two unique, locally valid perturbative approximations are complete. If 
desired, they can be combined to make a single uniformly valid approximation. 

While straightforward in theory, asymptotic matching is difficult in practice, particularly 
for an equation like the Navier-Stokes equation. However, it is still far simpler than alter- 
natives, such as iteratively solving the Oseen equations. Van Dyke's book is an excellent 
presentation of the many subtleties w hich arise in applying matched asymptotics to problems 
in fluid mechanics ( iVan Dykd . Il975l ). We now present the matched asymptotic solutions for 
Eqns. |9] and [12] These solutions result in the "state of the art" drag coefficients for both 
the sphere and the cylinder. 

a. Sphere 

Although Lagerstrom and Cole initially applied matched asymptotics to the problem of 
the sphere, the seminal work came in an elegant 1957 pa per by Proudman and Pearson 



( iLagerstrom and Cole 



1955 



Proudman and Pearson 



1957) . Ch e ster a nd Breach extended 



this paper via a difficult calculation in 1969 (IChester and Breach! . Il969l ). We summarize the 



results of both papers here. These workers used a perturbative solution in the Stokes regime 
of the form: 



V>(r, n) = V> + R&1 + R 2 log R1P2L + R 2 ^2 + R 3 log R + R 3 *p 3 + O (i? 3 ) 



(53) 



This rather peculiar perturbative form cannot be determined a priori. Rather, it arose in a 
fastidious incremental fashion, calculating one term at a time. The procedure of asymptotic 
matching required including terms like R 2 \ogR in the expansion; otherwise, no matching is 
possible. Note that matched asymptotics gives no explanation for the origin of these singular 
terms. 



is 



At this point, there are two unknown coefficients in the Oseen expansion, and one in the Stokes expansion. 
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The first step to finding a perturbative solution in the Oseen region is to define the Oseen 
variables: 

p = Rr, V(p,fi) = R 2 ^{r,(i) 

Part of the reason for this transformation can be understood via the derivation of the Oseen 
equation. The region where inertial effects become important has been shown to be where 
R\r\ ~ 0(1). Intuitively, the variable p = Rr is a natural choice to analyze this regime, 
as it will be of O (1). The precise reasons for the selection of these variables is a technique 
from boundary layer theory kn own as a dominant balance argument, which we will revisit 
later ([Bender and Orzagl . Il999l ) . 

The perturbative expansion in the Oseen region takes the form: 



V(p, n) = V + RVi + R 2 ^2 + R 3 log R^ 3 L + O (R 3 ) 



(54) 



Note that there is no R 2 log R term in this expansion; none is required to match with the 
Stokes expansion. As with the Stokes expansion, this form cannot be determined a priori. 

Proudman and Pearson completely solved for the Stokes' expansion through 0(R log-R), 
and partially solved for the O {R 2 ) term. They determined the Oseen expansion through 
O (R). Chester and Breach extended these results up to a partial solution for O (R 3 ) in the 
Stokes' expansion, and to O (R 3 \ogR) in the Oseen expansion. 

The exact form of these expansions is given in Chester and Breach. 16 Some aspects of 
these results have been seen before: The leading order in the Stokes' expansion (i/jq) is simply 
the Stokes solution (Eqn. [281) . In the Oseen expansion, ty is simply the formula for the 
uniform stream expressed in Oseen variables. The second term, is the rotational part of 
Oseen's solution (Eqn. [381) . 



Note that the expression for ip2 i n Proudman, is incorrect ( Proudman and Pearsonl Il957t ) . There is also a 



mistake in Chester and Breach ([Chester and Breach! I1969I ). Eqn. 3.5; the coefficient of cs should be 
not r~ 2 . 
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FIG. 10 (Color online) Drag on a sphere, experiment vs. matched asymptotic theory. Experimental 
and numerical results are plotted as in Figure [8] 

Both sets of authors then used their result for the Stokes expansion to calculate Co, 
which is given in Eqn. [SSJ 



C 



D 



"Stokes" "Oseen" 

40 V 3 & 



Proudman 

323N + ^R 3 log R+O (i? 3 ) 
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(55) 



Chester and Breach 

Here 7 is Euler's constant. This formula reproduces and extends nearly all earlier work. 



Eqn. [55] shows both the original resu 



contr ibutions of Chester and Breach (IChester and Breach 



ts of Proudman and Pear s on and the higher order 



1969 : 



Proudman and Pearson 



19571 ). The "Stokes" term is Stokes' original result (Eqn. [29]), which was rigorously justified 
by Oseen. The "Oseen" term is generally credited to Oseen (Eqn. 1391) . although it is really 
beyond the accuracy of his work, and is only justified by this calculation. 17 

Figure [10] compares the results of matched asymptotics (Eqn. [55]) with experimental 
data, numerical results, and the basic prediction of Oseen's equation (Eqn. [39]). This plot 
has been the source of some confusion. Maxworthy examined his data and concluded that 
Cd as computed b y Oseen and Gold stein (Eqn. [48]) were as good as any matched asymp- 



totics predictions (jMaxworthy 



196~1 . 



The calculations of Dennis and Le Clair, however, 



17 



Illingworth's unpublished result also justifies this term. 
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refute that claim, and demonstrate the systematic improvement that results from matched 
asymptotics. 

Neither is it immediately clear that the extra terms in Eqn. EH] due to Chester and 
Breach are actually any improvement on the work of Proudman and Pearson. Van Dyke 
notes, "This result is disappointing, because compar ison with experim ent suggests that the 
range of applicability has scarcely been increased." (IVan Dykd . Il975l ). and Chester himself 
remarks that "there is little point in continuing the expansion further." At very low Reynolds 
number, however, the results of Dennis "indicate that the expression of Chester and Breach 
gives a better appr oximation to the drag coef ficient than any other asymptotic solution until 



about [R = 0.3]." (IDennis and Walker 



19711 ) Figure [TU] shows the excellent low R agreement 



between Dennis' numerical results and Eqn. [531 

The prediction of matched asymptotics (Eqn. 155]) is close to Illingworth's second order 
Oseen theory (Eqn. 1521) . Close examination shows that the matched asymptotics results 
are slightly closer to the Dennis' calculations in the limit of low Reynolds number. Strictly 
speaking, these two theories should only be compared as r — > 0, and in this regime matched 
asymptotics is superior. This is not surprising, as the best matched asymptotic calculation 
is a higher order approximation than that of Illingworth. 

b. Cylinder 

In 1957, Kaplun applied matched asy mptotics to th e problem of the cylinder, and 
produced the first new result for Cd (jKaplunl . 119571 ). Additional stream function 
calculations (but without a drag coefficient) were done by Proudman and Pearson 
(jProudman and Pearson! . 119571 ). Kaplun's calculations were extended to higher order by 
Skinn er, whose work also explored the structure of the asymptotic expansions (jSkinnerl . 
19751 ). We summarize results from all three papers here. 

Near the surface of the cylinder, the Stokes' expansion applies, and the perturbative 
solution takes the following form. 



M, 9) = Mr, 0, 5) + Rip 2 {r, 6, 5) + R 2 ip 3 { r , 0,5) + (i? 3 ) 



(56) 



Here, 5 = S(R) is defined as in Eqn. HU What is remarkable about the structure of 
this expansion is its dependence on 5. To be precise, each function ip n is actually another 
perturbative expansion, in 5: 



Mr, 0, 5) = 5F nA (r, 0) + 5 2 F n , 2 (r, 9) + O (<5 3 ) 



(57) 
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This formulation is equivalent to an asymptotic expansion in terms of log R 1 , which is used 
by Proudman and Pearson: 

if fli m jkM) , F nt2 (r,e) ( 1 \ 

(log/?) 1 (logi?) 2 V(logi2) 3 y 
This form is much less efficient than that given in Eqn. [53, in the sense that more terms 
in the Stokes and Oseen expansions are needed to obtain a given number of terms in Co- 
For that reason, expansions in 5 are used here. 

This curious asymptotic form is necessitated by matching requirements. It is also the 
source of a number of bizarre complications. The first implication is that all terms in 
Eqn. [56] of O (R) and higher will be transcendentally smaller than any of the terms in the 
expansion for ipQ. This is true asymptotically, as R — » 0. The reason for this is that inertial 
terms never enter into any of the governing equations for the Stokes' expansion; they enter 
only through the matching process with the Oseen expansion. 

As with the sphere, the first step to finding a perturbative solution in the Oseen region 
is to transform into the relevant Oseen variables. In this case, 

p = Rr, ^(p,p) = /^(r,p) (59) 

The perturbative expansion which can solve the problem in the Oseen region has the same 
generic form as Eqn. [56j 

(p, 9) = %(p, 0, 5) + RViip, 8,5) + (R 2 ) (60) 

The functions ty n (p,6,5) can also be expressed as a series in S(R). However, the formula 
cannot be written down as conveniently as it could in Eqn. [53 The first two terms take the 
forms given in Eqn. [6TJ 

%(p, 9, 5) = F 0i0 (p, 0) + 6F 0A (p, 6) + <D (5 2 ) (61a) 
^(p, 9, 5) = 5- 1 F lt _ 1 (p, 9) + F 0i0 (p, 9) + O (5) (61b) 



Kaplun and Proudman both considered only terms of O (R°) in the Stokes' expansion. As 
R — > 0, this is an excellent approximation, as all higher terms are transcendentally smaller. 
In this limit, the Stokes expansion takes a particularly simple form: 



ip(r, 9) = ifti(r, 9, 5) = a n 5 n ilr logr — r H — J sin 9 

n=l ^ r ' 



(62) 
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Kaplun obtained terms up to and including n = 3. Proudman et al. also obtained expressions 
for the Oseen expansion, albeit expressed as a series in logi? -1 . Skinne r extended Kap lun's 



Stokes expansion to include terms up to O (5 3 ), O (R5), and O (R 2 S) (jSkinnerl . Il975l ). He 
obtained approximate solutions for the Oseen expansion, including terms up to O (5) and 
O (R). The lowest order solutions in the Oseen expansion are related to the expression for 
a uniform stream and the solution of Lamb (Eqn. |4"T1) . 

Using his solution, Kaplun computed a new result for the drag coefficient (Eqn. [631) which 
agrees with Lamb's result (Eqn. |4*2]) at lowest order. 

4-7T 

C D = — {5 - k5 3 ) (63) 

ri 

Here, k = J °° K (x)Ki(x) {x~ x I\{2x) — 4Ki(x)Ii(x) + 1) dx ~ 0.87. Skinner extended these 
results, showed that terms of O (R°) do not contribute to the drag, and calculated the first 
transcendentally smaller contribution, which is of O {R 1 )- His result is given in Eqn. [6H 

Cn = f (S - H 3 + O (* 4 ) -^(l- S - + (5 2 )) + O (i? 4 )) (64) 

The value of these new terms is questionable, and Skinner himself noted that they are 
likely negligible in comparison to the neglected terms of (9(5 4 ). Asymptotically this is 
unequivocally true. 

Figure [TT] compares the predictions of matched asymptotic theory with Lamb's result 
(Eqn. |4*2|) based on Oseen's equation. Although both theories agree as r — > 0, matched 
asymptotic results seem no better than Lamb's solution. The comparison is further compli- 
cated by the scatter in the different experiments; matched asymptotics agree more closely 
with Tritton's measurements, while Lamb's solution agrees better with Jayaweera's. We 
draw two conclusions from Figure [HJ Both theories break down for R > 0.1 and nei- 
ther theory is demonstrably more accurate. Even more disappointingly, Skinner's result is 
nowhere better than Kaplun's — it is actually worse at higher Reynolds numbers. 

Part of the problem with the matched asymptotics approach arises from the need for two 
expansions, in 5 and R. Because infinitely many orders of 5 are needed before any higher 
orders in R are relevant means that infinitely many terms in the Oseen expansion must be 
calculated before the second order term in the Stokes expansion. This is inefficient, and is 
the reason for Skinner's lack of success. 



A recent paper by Keller et al. solved this problem numerically ([Keller and Ward! Il996[ h 



They developed a method to sum all of the orders of 5 for the first two orders of R. Their 
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"beyond all orders" numerical results prove the importance of these higher order terms. 
When such terms are accounted for, the resulting Co is vastly improved from Kaplun's, and 
is superior to any of the analytic solutions discussed here. Interestingly, it seems to agree 
very well with the experiments of Tritton, although it is difficult to tell from the plot in 
their paper, which does not remove the leading order divergence. 



4- Other theories 



Amongst the community interested in settling velocities and sedimentation, there are 
many theoretical models of the drag on a sphere. These workers specify Cc as a function of 
R by means of a "sphere d rag correlation." An overview of these formula is given by Brown 



( IBrown and Lawlerl . 



20031 ) . These results are generally semi-empirical, relying on a blend 



of theoretical calculations and phenomenologically fit parameters to predict Co over a large 
range of Reynolds number. While practically useful, these results are not specific to low 
Reynolds numbers, and cannot be derived from the Navier-Stokes equations. They address 
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a different problem, and will not be further considered here. 



One other semi-empirical theory is due to Carrier ((Carrier 



19531 ). He argued that the 



inertial corrections in the Oseen equation were over-weighted, and multiplied them by a 
coefficient which he constrained to be between and 1. Consequently, his theory is in 
some sense "in between" that of Stokes and that of Oseen. He ultimately determined this 
coefficient empirically. 



5. Terminology 



Confusing terminology, particularly in the matched asymptotics literature, riddles the 
history of these problems. We previously detailed discrepancies in the definition of Co- In 
this section we explain the sometimes conflicting terms used in the matched asymptotics 
literature, introduce a convention which eliminates confusion, and also explain how some 
authors adopt different definitions of the Reynolds' number. 

Matched asymptotics literature discusses numerous perturbative expansions, each of 
which are valid in a different regime, or "domain of validity." Different authors use dif- 
ferent labels for these expansions. Most workers define the "inner" e xpans ion to be the 
expansion which is valid inside the boundary layer (IBender and Orzagi . Il999l ). A boundary 
layer is a region of rapid variation in the solution. The "o uter" expansion is valid ou tside of 
the boundary layer, where the solution is slowly varying (jBender and Orzagi . Il999l ). Prob- 
lems with multiple boundary layers require additional terminology. The outer expansion is 
based "upon the primary reference quantities in the problem," and the inner expansion is 
usually obtained by s tretching the ori ginal variables by dimensionless functions of the per- 
turbation parameter ( Van Dykel . 1975 ). The appropriate stretching, or scaling functions are 
obtained through a dominant balance analysis, which can be difficult. After this rescaling, 
the argument of the inner expansion will be of 0(1) inside the boundary layer. Accompany- 
ing these inner and outer expansions are "inner variables" , "outer variables" , "inner limits" , 
and "outer limits" . 

The low Reynolds number flow problems are complicated by the fact that so me authors , 



inclu ding Van Dyke, also define expansions on the basis of their physical location (jVan Dykd . 



19751 ). The "outer" limit is valid far from the solid body (|r| is large), and the "inner" limit 
is valid near the surface of the body (|r| ~ 1). 
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This is consistent with yet another definition, based on proximity to the origin of the 
chosen coordinate system. In a review paper Lagerstrom and Casten define the "inner 
limit" as bei ng "valid near the origin," an d the "outer limit" as being "valid except near 



the origin." (ILagerstrom and Casten 



1972) Part of their motivation for this new definition 



was to distinguish between the domain of validity of an expansion, and the limit process by 
which it is obtained. 

Finally, Kaplun refers to the inner and outer limits based on their correspondence to 
high Reynolds number flow (jKaplunl . 119571 ) . He identifies the Stokes' approximation as the 
"inner" limit, and Oseen's equation as the "outer" limit. 

Part of the confusion arises because of disagreements over the location of the boundary 
layer. Van Dyke claims that "it is the neighborhood of the point at infinity" , while Kaplun 
argues that the boundary layer is near the surface. Definitions referenced to the boundary 
layer disagree when there are disagreements about its location. 



To elimi nate this confusion, a pr e 
quent work ljKaplun and Lagerstrom 



erab l e alternative notation has eme rged from subse 



1957 



Proudman and Pearson 



1953). We follow this 



notation, defining the"Oseen" and "Stokes" expansions, which were used in the previous- 
section. The Oseen expansion is valid far from the surface, and is expressed in stretched 
coordinates. The Stokes limit is valid near the surface of the sphere, where r is small, and 
is expressed in the original variables. 18 

Matched asymptotics worke rs also discuss uniform approximations, intermediate expan 



sions, or composite expansions (IBender and Orzag 



1999 



Lagerstrom et al. 



1967 



Van Dykd . 



19751 ). The basic idea is that the Stokes and Oseen expansions can be blended together to 



form a single expression which is valid everywhere. This result reduces to the two original 
expansions when expanded asymptotically in the two limits. How to calculate a uniform 
expansion is discused below. 

There are also minor differences in the definition of the Reynolds number, R. Some 
authors define R based on the diameter of the solid, while others base it on the radius. This 
factor of 2 can be difficult to track. We define the Reynolds nu mber using the radius o f 
the fixed body: R = lu^a/v. It is worth noting that Kaplun (ILagerstrom et all 119671 ). 



18 Van Dyke's book is not consistent in relating "inner" and "outer" expansions to the Stokes and Oseen 
expansions. 
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Tomoti ka (ITomotika and Aoil . Il950f) . Goldstein ([Goldstein! I1965T 1. Liebster (jLiebsterl . 119271 ). 
Thorn (jThoml . Il933l ) and Tritton (ITrittonl . Il959l ) all use the diameter. 



C. Uniformly valid approximations 

As mentioned previously, the inner and outer expansions may be combined into a sin- 
gle, uniformly valid approximation, which is applicable everywh ere. For a function of on e 



variable, the uniform approximation is constructed as in Eqn. [65] (jBender and Orzag 



l/uniform (*^) 2/outer(*^) Z/inner (*^) 2/overlap (^) (65) 

y ovcr i ap (x) consists of the common "matching" terms between the inner and outer expansions. 
Kaplun demonstrates that I/uniform^) — ► 2/(^) as the expan sion variable -R — > 0, i.e . the 



uniform approximation tends to the exact solution everywhere (ILagerstrom et all , 119671 ). To 



be more precise, if the matched asymptotics solution is constructed to O (R 1 ), then 

lim y(x) - y miiiovm (x) ~ O (R 1 ) 

As a matter of practice, calculating the uniform solution is mechanistic. First, express 
the inner and outer expansions in the same coordinates; in our case, express the Oseen 
expansion in Stokes variables. 19 Alternatively, one can express the Stokes expansion in 
Oseen variables. Next, express both solutions as a power series in the expansion parameter, 
R. By construction the Stokes expansion is already in this form but the transformed Oseen 
expansion is not, and must be expanded to the same power in R as the Stokes solution. 

From these two power series we can identify the "overlap" function, ^overlap- This function 
consists of the terms which are in common between the two expansions, and is usually 
obtained by inspection. Of course, ^overlap is only valid to the same order as the original 
matched asymptotics solution, and higher order terms should be discarded. The uniformly 
valid approximation is then obtained using ^overlap and Eqn. 1651 



19 Note that this transformation affects both the radial coordinates and the stream function, and that it 
differs for the sphere and cylinder. 



57 



1. The correct way to calculate Cd 



Proudman and Pearson argue that "uniformly valid approximations per se are not usually 
of much physical interest ... In the present problem, for instance, i t is the Stokes expansion 



that gives virtually all the physically interesting information." (IProudman and Pearson 



19571 ) All matched asymptotics calculations are based solely on the Stokes expansion, and 



are therefore influenced by the Oseen expansion only via the boundary conditions. For 
instance, the drag coefficient is calculated using only the Stokes' expansion. Other properties 
of the stream function, such as the size of the dead water wake directly behind the sphere 
or cylinder, are also calculated using the Stokes' expansion. 

In this section we argue that this approach is incorrect, and that uniformly valid approx- 
imation should be used to calculate all quantities of interest. By adopting this viewpoint, 
we obtain new results for Cd, and demonstrate that these drag coefficients systematically 
improve on previous matched asymptotics results. 

Matched asymptotics workers argue that the drag coefficient is calculated at the surface 
of the solid (Eqns. [2U [T9"j) . where r = 1. Since the Oseen solution applies for large r, the 
Stokes solution applies for small r, and the Stokes solution ought to be used to calculate 
Co- In fact, by construction, any uniformly valid approximation must reduce to the Stokes 
expansion in the limit as Rr — > 0. 



Curiously, proponents of the Oseen equation argue conversely (IFaxenl . Il927 : 



Happel and Brenner! . Il973l ). They claim that because the Oseen expansion happens to ap- 



ply everywhere, it should be used to calculate all sorts of quantities of interest, includ- 
i ng Cn- In fact, Hapel an d Brenner wrote a book essentially devoted to this premise 



(IHappel and Brenner! . 



19731 ). In fairness, it must be mentioned that all of these authors 



were well aware of their choices, and motivated their approach pragmatically: They ob- 
tained useful solutions to otherwise intractable problems. 

In reality, both approaches converge to the exact solution for suitably small Reynolds' 
numbers. However, for small but non-infinitesimal R, the best estimate of derivative quan- 
tities such as Cd is obtained not by using the Stokes expansion, but by using a uniformly 
valid approximation calculated with both the Stokes and Oseen expansions. Such a drag 
coefficient must agree with results derived from the Stokes expansion as Rr —>■ 0, and it 
can never be inferior. Moreover, this approach makes determination of the drag coefficient's 



58 



accuracy straightforward; it is determined solely by the accuracy of the uniform expansion, 
without any need to be concerned about its domain of applicability. 

We now calculate the drag coefficients for both the sphere and the cylinder using uni- 
formly valid approximations, using previously published inner and outer expansions. These 
corrections are small but methodologically noteworthy, and are absent from the existing 
literature. 

a. Cylinder 

Although the state-of-the-art ma tched asympto tics solutions are due to Kaplun, it is more 
convenient to use stream functions (IKaplunl . I1957I ). Skinner conveniently co mbines previou s 
work, providing a concise summary of Stokes and Oseen stream functions (jSkinnerl . Il975l ). 
We base our derivation of a uniformly valid approxim ation on the results in his paper. The 



19751 ). 



Stokes expansion is given by Eqn. [66] (ISkinner 

^(r, 9) = ~ (5 - kS 3 + O (<5 4 )) \2r log 
The Oseen Expansion is given by Eqn. [671 



r — r + 



ij sinO + 0(R 1 ) (66) 



/ oo 

V(p,6) = Lsm9-5j2^n (0 ^smn9 + 0(5 2 ) + 0(R 1 ) 



(67) 



n=l 



With these results, creating the uniform approximation and calculating Co is straight- 
forward. The only subtlety is the sine series in Eqn. [67J However, Eqn. [21] tells us that, for 
the purposes of calculating the drag, only the coefficient of sin 9 matters. We calculate the 
overlap between the two functions by expanding Eqn. [67] about p = 0. The result is given 
by Eqn. ESI 

VWria P (r, 9) = 8~ (2 logr - 1) sin 9 + O {5 2 ) + O {R 1 ) (68) 

Combining this with the Oseen and Stokes expansions, we obtain the uniformly valid ap- 
proximation given by Eqn. [691 



V^uniform (f, $) 



r + 5 

oo 
n=2 



r0i( 



1 

2r 
Rr\ r 
~2~ J n 



rR, 



+ k5 3 



r 1 
r log r 

2 & 2r 



sin 9 



smn9 + 0(5 2 ) + O (R 1 ) 



(69) 



By substituting this result into Eqn. EH we obtain a new result for Cd- 

tx5 (24 - 32M 3 + QR 2 <f)'i{R/2) + i? 3 0' 1 "( J R/2)) 



C 



D 



8R 



(70) 
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FIG. 12 (Color onli ne) Drag on a cylinder 



asymptotics results (jJavaweera and Masonl . 



comparing a uniform ly valid calculations and matched 



1965 



Trittonl . 



19591 ) 



Fig. [12] compares Eqn. [70] with Kaplun's usual result (Eqn. 1631) . The new drag coefficient 
(Eqn. [701) is a small but systematic improvement over the results of Kaplun. Because they 
are asymptotically identical up to O (<5 4 ) and O(R), they agree as R — > 0. However, at 
small but non-infinitesimal R, our new result is superior. Comparing Figures [12] and [TU we 
can also see a second surprise: The new result betters Skinner's Cp, even though they were 
based on the same stream functions. If Skinner had used a uniformly valid approximation, 
his result would not have misleadingly appeared inferior to Kaplun's. 

b. Sphere 

As with the cylinder, calculating Cd from a uniformly valid expansion yields an improved 
result. However, there is a substantial difference in this case. Although matched asymptotics 
calculations have been done through O (i? 3 ) in Eqn. [53] and O [R? log R) in Eqn. [5H the 
higher order terms in the Oseen expansion are impossible to express in a simple analytic 
form. Asymptotic expressions exist (and have been used for matching), but these cannot 
be used to construct a uniformly valid expansion. Consequently, we can only compute the 
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uniform expansion through O(R), and its predictions can only be meaningfully compared 
to the first two terms in Eqn. 1551 

The solutions for the Stokes and Oseen e xpansions are given in Chester and Breach, and 



are quoted here ([Chester and Breach! . Il969l ). The Stokes expansion: 



^(r,/i) = -^(2r 2 -3r+^Q 1 (^)-R^((2r 2 -3r + ^Q 1 {fx)- 

^2r 2 - 3r + 1 - i + Q 2 (jm) j + O (R 2 logR) (71) 
The Oseen expansion: 

¥(p, n) = - P 2 QM - R\ (1 + /i) (l - e^ 1 ^) + O {R 2 ) (72) 

By taking the p — > limit of Eqn. [721 we can calculate the overlap between these two 
expansions. The result is given in Eqn. [731 

T/Veria P (r, //) = £ (12 - 8r) + ^ (3rQ 2 (/i) - 3^^)) + (i? 2 ) (73) 

o o 
Eqns. [731 E21 and [7T] can be combined to form a uniformly valid approximation: 

VWorm(r, //) = Vfa /") - VWrlap(r, fl) + ^ ^ ^ + O (R 2 log fl) (74) 

Due to the e - "^ 1- **) term, we cannot use the simple expression for Cd (Eqn. [26]) . Instead, 
we must use the full set of Eqns. [HI [241 and [251 After completing this procedure, we obtain 
a new result for Cd, given by Eqn. [75j 

g / -2R 

C D = — t ( 40e fl (1728 + 1140.R + 335# 2 + 56i? 3 + 6R 4 ) - 60R (1 + R) 

R \320R 3 V v ; v ; 

+e 2i? ( - 69120 + 23580i? - 2420# 2 + 20(10 + tt)R 3 + 10(18 - 7r)i? 4 - 8R 5 



e' R / 2 7r 



h(R/2) \ 



- 3R6 ))- £ (75) 

This result is plotted in Figure [131 Asymptotically, it agrees with the matched asymp- 
totics predictions to 0(1), as it must, and reproduces the 3/8R "Oseen" term. As R 
increases, however, the uniform calculation becomes superior to the first two terms of the 
matched asymptotic Cd- Although it is a much higher order solution than either of the 
other two results, we show the full matched asymptotics prediction for comparison. 
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III. THE RENORMALIZATION GROUP APPLIED TO LOW R FLOW 
A. Introduction to the renormalization group 

In 1961, Lagerstrom proposed the first of a number of "model problems", ordinary differ- 
ential equations which exhibited many of the same asymptotic features as the low Reynolds 
number problems. They were used to study and develop the theory of matched asymptotic 
expansions. The mathematical solution of these problems is closely analogous to the actual 
solutions of the Navier-Stokes equations. 



A re view of these equations, and o 



their matched asymptotic solutions, is given in Lager- 



strom (ILagerstrom and Castenl . 119721 ). The relevant models can be summarized by the fol- 



lowing equation: 

d 2 u n — ldu du . f du\ 2 „ 

+ — +U— + 81 — ) =0 (76) 



dx 2 x dx dx \dx / 
This ODE is subject to the boundary conditions u(e) = 1, u(oo) = 0. In this equation, 
n corresponds to the number of spatial dimensions (n = 2 for the cylinder, n = 3 for the 
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sphere). 5 = characterizes incompressible flow, and 5=1 corresponds to compressible flow. 
This equation is similar to the Navier-Stokes equations expressed in Oseen variables. There 
are fundamental differences between the structure of the incompressible and compressible 
flow equations. 

These model pr oblems are p osed in Hinch, albeit in terms of "Stokes" (rather than 



'Oseen") variables ( Hinch! . 



1991]) . Hinch begins by examining the model describing incom- 



pressible flow past a sphere. He next examines incompressible flow past a cylinder, which he 
calls "A worse problem." Finally, he treats compressible flow past cylinder, which he dubs 
"A terrible problem." 

These problems, which have historically been the proving ground of matched asymptotics, 
were recently solved using new R e norm alization Group (RG) techniques in two papers by 



Chen et al. (IChen et al. 



1994b 



19961 ). These techniques afford both quantitative and 



methodological advantages over traditional matched asymptotics. The RG approach derives 
all of the subtle terms (e.g., R 2 \ogR) which arise during asymptotic matching, demonstrat- 
ing that origin of these terms lies in the need to correct flaws inherent in the underlying 
expansions. Moreover, RG does not require multiple rescalings of variables, and its results, 
while asymptotically equivalent to those of matched asymptotics, apply over a much larger 
range (e.g., they extend to higher R). 

In particular, Chen et al. solved Hinch's first model, which describes incompressible flow 
past a sphere (n = 3, 5 = 0), a s well as the model for both kinds of flow past a cylinder 



n 



2, 5 = 0,1) (jChen et al 



1994b 



19961 ). In a notation consistent with Hinch, they 



termed these models the "Stokes-Oseen caricature" and the "terrible problem." 

The dramatic success of the RG techniques in solving the model problems inspired their 
application to the original low Reynolds number flow problems. That is our primary purpose 
here, as the low Reynolds number problems are the traditional proving ground for new 
methodologies. We will show that the RG techniques perform well when applied to these 
problems. RG produces results superior to and encompassing the predictions of matched 
asymptotics. More importantly, the RG calculations are considerably simpler than matched 
asymptotics, requiring half the work. 

The utility of the RG approach is most easily seen through an example, which will also 
provide a framework for understanding the analysis presented in su bsequent sections. Several 



pedagogical examples can also be found in the references (e.g., (IChen et all . Il994bl . 



1996: 
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Goldenfeldl . Il992l ; lOond . |2000| )). We begin here with an analysis of the most complicated 
model problem, the "terrible problem," which caricatures compressible flow past a cylinder. 



1. Detailed analysis of the "terrible problem" 

Although the "terrible problem," is solved in a paper by Chen et al., we re-examine it here 
in considerably more detail, as its solution is closely analogous to those of the low Reynolds 
number flow problems. This switchback problem is exceptionally delicate 20 , requiring the 
calculation of an infinite number of terms for the leading order asymptotic matching. 

There are pitfalls and ambiguities in applying RG techniques, even to the "terrible prob- 
lem," which while terrible, is considerably simpler than the real low Reynolds number prob- 
lems. Understanding these subtleties in this simpler context provides essential guidance 
when attacking the Navier-Stokes' equations. 

We want to solve the ODE given in Eqn. I77al subject to the boundary conditions I77bl 

This equation can be derived from Eqn. [76] by setting n = 2, 5 — 1, and transforming to the 

"Stokes" variables, r = x/e. Unlike Eqn. [76], Eqn. [77] is obviously a singular perturbation 

in e, which has been removed from the boundary conditions. The last term in the equation 

vanishes when e = 0. 

d 2 u{r) ^ \du{r) ^ / du(r)\ 2 ^ du{r) ^ 
dr 2 r dr \ dr J dr 

u(l) = 0, u(r = oo) = 1 (77b) 

This problem cannot be solved exactly, although numerical solution is straightforward. 
Trouble arises due to the boundary layer 21 located near r = oo. RG analysis requires that 
we work in the "inner" variable for our approximation to capture the correct behavior near 
the boundary layer 22 . This requirement may also be qualitatively motivated by arguing 
that one must choose coordinates to "stretch out" the boundary layer so that it can be well 
characterized by our approximate solution. 

To determine the appropri ate change of variables, w e need to analyze Eqn. ( 1771) using a 



dominant balance argument (iBender and Orzagj . Il999l ). As it stands, the first three terms 



20 Hinch notes, "It is unusual to find such a difficult problem ..." (jHinchl . Il99ll ). 

21 A boundary layer is a region of rapid var iation in the solution, y(t) . 

22 Here we use "inner" in the usual sense ([Bender and Orzad . Il999l ). For further discussion, see Section 
III.B.5I 
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of Eqn. (I77ap will dominate, since e is small. The rescaling x = er yields inner Eqn. 
This, of course, is the same equation originally given by Lagerstrom (Eqn. 176]) . 

d 2 u{x) ^ 1 du(x) ^ / du(x)\ 2 ^ du{x) ^ 
dx 2 x dx \ dx J dx 

u{e) = 0, u{x = oo) = 1 (78b) 



The next step in the RG solution is to begin with the ansatz that the solution to Eqn. 
rSj) can be obtained from a perturbation expansion (Eqn. [791) . We fully expect this ansatz 
to fail, since we have a singular perturbation in our ODE. We therefore refer to this starting 
point as the naive perturbation expansion. 

u(x) = u (x) + eux{x) + e 2 u 2 {x) + C(e 3 ) (79) 

Collecting powers of e, we obtain differential equations for u (x), Ui(x), etc: 

O (e°) : + u (x)u' (x) + u' (x) 2 + u0"(x) = (80) 



O (e 1 ) : M!M + — + m Mi + 2m m 1 + uj = (81) 
O (e 1 : ^2^0 + ■u 1 Mi + M0M2 + + 2m m 2 H h m 2 = (82) 



x 

u 2 

X 

a. O(e ) solution 

The first complication of the terrible problem arises when we attempt to solve Eqn. 
[80], a nonlinear ODE. Although one solution — Uq(x) = A — is seen by inspection, an 
additional integration constant is not forthcoming, and our solution to the O(e ) problem 
cannot satisfy both of the boundary conditions (Eqn. I78bl) . The resolution to this quandary 
is simple: Ignore the problem and it will go away; continue constructing the naive solution 
as if Uq(x) = A Q were wholly satisfactory. The qualitative idea is that the O(e ) solution is 
the uniform field which we have far from any disturbance source. Why is this acceptable? 

The RG method is robust against shortcomings in the naive expansion. We know that 
singular perturbation problems cannot be solved by a single perturbation expansion. We 
therefore expect problems, such as secular behavior, to arise in our solution for the naive 
expansion. RG techniques can be used to remove t hese flaws from the perturbative solution, 



turning it into a uniformly valid approximation (jChen et all Il996l ). It does not matter 



whether these defects arise from an incomplete solution for Uq,(x), the intrinsic structure of 
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the equation, or a combination of the two. To solve the terrible problem (and later the low 
Reynolds number problems), we must exploit this flexibility. 

For subsequent calculations, there are two ways to proceed. First, we may retain A as an 
arbitrary constant, one which will ultimately be renormalized in the process of calculating a 
uniformly valid approximation. Alternatively, we may set A = 1, satisfying the boundary 
condition at x = oo. 23 This unconventional approach to the RG calculation effectively 
shifts the freedom that usually comes with the O(e ) constants of integration into the 0(e l ) 
solution. This artifice greatly simplifies subsequent calculations, and is invaluable in treating 
the Navier-Stokes equations. Moreover, these two approaches are equivalent, as we now 
show. 

b. 0{e 1 ) solution 

If u (x) = A , Eqn. [811 simplifies to Eqn. [83j 

d 2 Ui ( 1 ^ \ du\ . . 

dx 2 \x / dx 

The solution is: Ui(x) = B + Biei(A x), where e n (x) = e~H~ n (it. Notice that the first 
term is linearly dependant on the Un(x) solution. There are ma ny opinions regarding how to 



utilize this degree of freedom (jKunihird . Il995l ; IWoodrufil . Il995l ) . In our approach, one is free 
to choose the homogeneous solutions of Uo, U\, etc. for convenience. The only constraint 24 is 
that the "naive" solution (Eqn. [79]) must have a sufficient number of integration constants 
to meet the boundary conditions. In this example, that means two constants of integration. 

Different choices of particular solutions will ultimately result in different approximate 
solutions to the ODE. However, all of these solutions will agree within the accuracy limita- 
tions of the original approximation (in this case the naive expansion). This can be shown 
explicitly. In this example, as in the low Reynolds number problems, we choose a particular 
solution which simplifies subsequent calculations. Setting B = (note that this is not the 
same as a redefinition of the constant A ), we obtain the solution: 

u(x) = A + eBxexjApx) +0(e 2 ) (84) 

divergent as x— *0 

The second term in Eqn. [841 diverges logarithmically as x — > 0. One may argue that this 
divergence is irrelevant, since the range of the original variable is r 6 [1, oo), and numerical 



23 Meeting the boundary condition at x = e results only in the trivial solution uq(x) = 0. 

24 Of course the solution must also satisfy the governing equation. 
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solutions demonstrate that the solutions to Eqn. [TTJin [1, oo) diverge when extended to r < 1. 
But the argument that the divergence in Eqn. [HH is an intrinsic part of the solution (and 
therefore should not be considered problematic) is incorrect. Although the original variable, 
r, is limited to r G [1, oo), the transformed variable, x = er, has the range x G [0, oo). This 
occurs because there are no restrictions on the lower limit of e. The divergence exhibited 
by the second term of Eqn. [HH must be removed via renormalization in order to turn the 
flawed naive solution into a uniformly valid approximation. 

This divergence arises for two reasons. First, we are perturbing about an O(e ) solution 
which is deficient; it is missing the second integration constant (and concomitant funda- 
mental solution). More fundamentally, Eqn [7JJ] attempts to solve a singular perturbation 
problem with a regular expansion, an approach which must fail. The RG technique solves 
these problems by restructuring the naive expansion and eliminating the flaws in uq(x). 

Although A is simply a constant of integration when e = 0, it must be modified when 
e 7^ 0. We will absorb the divergences into a modification, or renormalization, of the constant 
of integration Aq. Formally, one begins by "splitting" the secular terms, replacing ei(Aox) 
by ei(A x) — ei(A r) + ei(A T), where r is an arbitrary position. This results in Eqn. ESJ 

u{x) = A + eB 1 (e 1 (A x) - e 1 {A r) + e 1 (A r)) + 0{e 2 ) (85) 

Since r is arbitrary, it can be chosen such that ci(Aqx) — ei(A$T) is non-secular (for a given 
x). The divergence is now contained in the last term of Eqn. ([85]) . and is exhibited as a 
function of r. 

It is dealt with by introducing a multiplicative renormalization constant, Z\ — 1 + 
ai(r)e\ and then renormalizing A as A = ZiA q (t). 25 The coefficients Oj(r) can then 
be chosen. 26 order by order so as to eliminate the secular term in Eqn. f[85]) . Substituting, 
and choosing a± to eliminate the final term of Eqn. [85] we obtain 



u(x) = A (t) + eB 1 (e 1 (A (r)x) - e 1 (A (r)r)) + 0(e 2 ) (86) 

MrT~ (87) 

25 Aq is the only constant which can be renormalized to remove the divergences, as B\ is proportional to 
the secular terms. 



Where a\ satisfies 

aM = AM 



2i, 



Note that the coefficients must also be independent of x. 
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Note that to obtain Eqn. [86] we needed to expand e± about e = 0. Unusually in this 
equation, the renormalized constant (Ao(t)) appears in the argument of the exponential 
integral; this complicates the calculation. We will later show how to avoid this problem by 
restructuring our calculations. 

Qualitatively, the idea underlying Eqn. [HH] is that boundary conditions far away (from 
x = e) are unknown to our solution at i > c, so that Aq is undetermined at x = r. RG 
determines Aq in this regime through the renormalization constant Z\ (which depends on r). 
Afterward there will be new constants which can be used to meet the boundary conditions. 

The RG condition states that the solution u(x) cannot depend on the arbitrary position 
t. This requirement can be implemented in one of two ways. First, since d T u(x) = 0, apply 
d T to the RHS of Eqn. [HS and set the result equal to zero: 



A' (t) + eB 1 ( 6 -^-^ + 4H4 (e- Ao(T)T ~ e~ A ^)) + 0(e 2 ) = 



(88) 



The next step in RG is to realize Eqn. [HS implies that A' (t) ~ O(e). Retaining only terms 
of 0(e), we obtain: 

^i +cBl (L^iy o ^= (89) 

In principle, we simply solve Eqn. [89] for Aq(t). Unfortunately, that is not possible, 
due to the presence of Aq(t) in the exponential. This complication also occurs in other 
switchback problems, as well as in the low Reynolds number problems. Eqn. [89] can be 
solved by an iterative approach: Initially set e = 0, and solve for A (t) = a , a constant. 
Next substitute this result into the 0(e) term in Eqn. [89] solving for A (r) again: 

A {t) = a + eB iei {a T) (90) 

In this solution, we have a new integration constant, o.q. Having obtained this result, we 
again must exploit the arbitrary nature of r. Setting t = x, and substituting into Eqn. [86], 
we obtain: 

u (x) = a + eBiei(a x) + 0(e 2 ) (91) 

But this is identical to the original solution (Eqn. 183]) ! What have we accomplished? 
This renormalized result is guaranteed to be a uniformly valid result, for Vx. The renor- 
malization procedure ensures that the logarithmic divergence in Eqn. [91] is required by the 
solution, and is not an artifact of our approximations. Obtaining the same answer is a 
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consequence of solving Eqn. [88] iteratively. Had we been able to solve that equation exactly, 
this disconcerting coincidence would have been avoided. 

We obtain the final solution to Eqn. I77al by applying the boundary conditions (Eqn. 
I78bl) to Eqn. [HU a — 1, Bi — — l/(eei(e)). Lastly, we undo the initial change of variables 
(r = x/e), yielding the result given in Eq n. |9"2"1 As shown in Chen et al., this is an excellent 
approximate solution (IChen et all Il996l ). 

«(r) = 1-^ + (e 2 ) (92) 

Furthermore, if we expand the coefficient B\ = — l/(eei(e)) for e — > + , Bi(e)/e ~ 
— l/ln(l/e) — 7/ln 2 (l/e). These logarithmic functions of e are exactly those which are 
required by asymptotic matching! These "unexpected" orders in e make the solution of 
this problem via asymptotic matching very difficult. They must be deduced and introduced 
order by order, so as to make matching possible. In the RG solution, they are seen to arise 
naturally as a consequence of the term l/ei(e). 

There are several other equivalent ways to structure this calculation. It is worthwhile to 
examine these (and to demonstrate their equivalence), in order to streamline our approach 
for the low Reynolds number problems. 

The first variation occurs in how we apply the RG condition. Rather than applying d T 
to Eqn. [86], we may also realize that the original constants of integration, Aq = Zi(t)Aq(t), 
must be independent of r. Hence the "alternative" RG equation: 

dA _ d(Z 1 (T)A Q (r)) _ q 
dr dr 

Substituting Z 1 = 1 + e (- J Bie 1 (rA (r)(l + J™1 a i( T ) el ))) / A o( T ) + (e 2 )> one obtains: 

A' (r) + eB 1 ( + ^e~ A ^A + O (e 2 ) = (93) 

Because this implies A' (t) ~ O (e 1 ), Eqn. [931 simplifies to Eqn. EH] (to within O (e 2 )), and 
these two methods of implementing the RG condition are equivalent. 

In addition to this dichotomous implementation of the RG condition, there is yet another 
way to structure the analysis from the outset: We set A = 1 in the zeroth order solution, 
and rely on the robustness of the RG approach to variations in our perturbative solution. 
With this uq(x) solution, there is no longer any freedom in our choice of u\(x) integration 
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constants — both are needed to meet boundary conditions. In this approach, our naive 
perturbative solution is: 

u{x) = 1 + e(B + £iei(aO) + O (e 2 ) (94) 

divergent 

Proceeding as before, replace e\{x) by e\{x) — ei(r) + ei(r): 

u{x) = l + e(B + B 1 ( ei (x) - ei(r) + ei(r))) + O (e 2 ) 

Again introduce renormalization constants {Z\ = 1 + Y^Li a «( r ) e \ Z 2 = 1 + Y^Li bi{j)e l ), 
and renormalize Bq, B\ as Bq = ZiBq(t) and B\ = ZiB\{j). In fact, only Bq needs 
to be renormalized, as the Bi term multiplies the secular term and consequently cannot 
absorb that divergence. This can be seen systematically by attempting to renormalize both 
variables. With an appropriate choice of coefficients, ai = — Si(r)ex(r) and b\ = 0, the final 
term in the last equation is eliminated. b± = demonstrates that B\ does not need to be 
renormalized at O (e 1 ). The resulting equation is given in Eqn. [951 

u(x) = 1 + e (B (t) + 5i(r) ( ei (x) - e x (r))) + (e 2 ) (95) 

We did not actually need to determine ai or bi in order to write the above equation; it 
could have been done by inspection. Determination of these quantities is useful for two rea- 
sons. First, it helps us see which secular terms are being renormalized by which integration 
constants. Secondly, it allows the second implementation of the RG condition which was 
described above. This can sometimes simplify calculations. 

Using the first implementation (requiring d T u(x) = 0), and using Eqn. [951 we obtain: 

B' (r) + B[{r) ( ei (x) - ei (r)) + B x (r)^ = O (e 1 ) (96) 

This can only be true Vx if B[(t) = 0, or -Bi(t) = /3 2 , a constant (as expected). Knowing 
this, we solve for B (t) — (3\ + ^2^1 (r). Substituting this result into Eqn. [951 an d setting 
t = x, we obtain the renormalized solution: 

u{x) = l + e{(3 1 +(3 2 e 1 {x)) (97) 

The boundary conditions in Eqn. 1 78b I are satisfied if /?i = and /3 2 = — l/(eej.(e)). Returning 
to the original variable (r = x/e), we obtain: 

u(r) = l- G j^- + (e 2 ) (98) 
ei(ej 
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This is identical to Eqn. [92J, demonstrating the equivalence of these calculations. The latter 
method is preferable, as it avoids the nonlinear RG equation (Eqn. [89]) . We will use this 
second approach for analyzing the low Reynolds number problems. 

The RG analysis has shown us that the logarithmic divergences present in Eqn. |HH are 
an essential component of the solution, Eqn. [53 However, we must work to O (e 2 ) in order 
to see the true utility of RG and to understand all of the nuances of its application. 

c. O (e 2 ) solution 

We base our treatment of the O (e 2 ) on the second analysis presented above. Through 
O (e 1 ), the naive solution is: uq{x) = 1, ui(x) = B + Bie\(x). Substituting into Eqn. [821 
we obtain the governing equation for u 2 (x): 



This has the same homogeneous solution as Ui(x), u 2 (x) = C + C\ei{x). A particular 
solution is: 



As discussed previously, we are free to choose Co, G\ to simplify subsequent calculations. 
The constants B Q , Bi are able to meet the boundary conditions, so there is no need to retain 
the O (e 2 ) constants: We choose Cq = 0, C\ = 0. In this case, the differing differing choices 
of C , Ci correspond to a redefinition of B , B\ plus a change of O (e 3 ), i.e. B = B + eC . 27 
Our naive solution through O (e 2 ) is thus: 



The underlined terms in this expression are divergent as x — > 0; the doubly underlined 
term is the most singular (~ ln(x) 2 ). RG can be used to address the divergences in Eqn. 
I10QL However, there is a great deal of flexibility in its implementation; while most tactics 
yield equivalent approximations, there are significant differences in complexity. We now 
explore all of the organizational possibilities in the terrible problem, an exercise which will 
subsequently guide us through the low Reynolds number calculations. 

27 This was not true at the previous order. 




(99) 





(100) 
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The first possibility is to treat only the most secular term at O (e 2 ). The doubly under- 
lined term dominates the divergent behavior, and contains the most important information 
needed for RG to construct a uniformly valid approximation. The approximation reached 
by this approach is necessarily inferior to those obtained utilizing additional terms. However 
it is nonetheless valid and useful, and eliminating most of the O (e 2 ) terms simplifies our 
calculations. 

Discarding all O (e 2 ) terms except the doubly underlined term, we begin the calculation 
in the usual manner, but come immediately to the next question: Ought we replace e\(x) by 
e\{x) — e 2 (r) +e 2 (r) or by (ei(x) — ei(r)) 2 + 2ei(x)ei(r) — e 2 (r)? Each option eliminates the 
divergence in x, replacing it with a divergence in r. Both merit consideration. Beginning 
with the latter, the renormalized perturbative solution is: 

'1 



u(x) = l + e (B (t) + B 1 (t) {e x {x) - ei (r))) - e 2 (^iOT ( e i(*) ~ e i( r )Y 

+e 2 (less divergent terms) + O (e 3 ) (101) 

Applying the RG condition (d T u{x) = 0) results in a lengthy differential equation in r. 
Because we want our solution to be independent of x, we group terms according to their x 
dependence. Recognizing that B[(t) ~ 0(e 1 ), B (t) ~ 0{e l ), and working to C(e 3 ), we 
obtain two equations which must be simultaneously satisfied: 

B[( T ) - ' e " T - B ' (r) = O (e 3 ) (102a) 



r 



(102b) 



Eqn. I102al has the solution 



Substituting this result into Eqn. Il02b| and solving, we obtain the result 



flo(r)^ft + ' 1 ' (a+ ; e ' (T)) + 0(^) 

Both (3q and j3\ are constants of integration which can be later used to meet the boundary 
conditions. Substituting these solutions into Eqn. I101[ setting r = x, disregarding terms of 
O (e 2 ) and higher we obtain the renormalized solution: 

u(x) = 1 + e I p H 1 + O (e ) (103) 
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Choosing /3 and /3i to satisfy Eqn. I78bl results in Eqn. 11041 

u(x) = In (e + (1 ~^ 6 e ) l(X) ) + ^ (e 2 ) (104) 
Expressing this in the original variable (r = x/e), results in the final answer (Eqn. 11051) . 

u(r) = In fe + ii^iM^ + f£\ (105 ) 
V ei(e) J 

This is the solution previously obtained by Chen et al., albeit with a typographical 



error corrected (IChen et all Il996l ). We will now revisit this analysis, using the alternative 
"splitting" of the most secular term in Eqn. 11001 but not yet considering less secular (or 
non-secular) terms of O (e 2 ). 

If we replace replace e 2 (x) in Eqn. llOOl by e\(x) — e 2 (r) + e 2 (r), we obtain the new naive 
expansion given by Eqn. 11061 

u(x) = 1 + 6 (Bb(r) + B 1 (r) {e x {x) - ei (r))) - e 2 (^i(r) 2 (e 2 (x) - e 2 (r))) 

+e 2 (less divergent terms) + O (e 3 ) (106) 

We now repeat the same calculations: 

1. Apply the RG condition (d T u(x) = 0). 

2. Group the resulting equation according to x dependence. This will result in two 
equations which must be satisfied independently. 

3. Discard terms of O (e 3 ), observing that B' (t), B[(t) must be of O (e 1 ). 

4. Solve these differential equations simultaneously for B (r), Bi(r). 

5. Substitute these solutions into the original equation (i.e. Eqn. 11061) . and set t = x. 

6. Choose the integration constants in this result to satisfy Eqn. I78b[ 

7. Obtain the final solution by returning to the original variable, r = x/e. 

For Eqn. 11061 steps 1-4 result in the following solutions for our renormalized constants: 
B 1 (t) = (3 X + O (e 2 ), B (t) = /3o + PMr) - e(3le\{T)/2 + O (e 2 ). Completing step 5, we 
obtain the renormalized result: 

u{x) = 1 + 6 (A, + PMx)) - e 2 ^ + O (e 2 ) (107) 
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This is identical to our starting point, Eqn. 11001 (retaining only the most secular terms). 
This should no longer be surprising, as we observed the same phenomena in the O (e 1 ) anal- 
ysis (Eqn. liTTj) . However, it is worth noticing that we obtained two different results (Eqns. 
11041 1107p depending on how we structured our RG calculation. This apparent difficulty is 
illusory, and the results are equivalent: Expanding Eqn. 11031 for small e reproduces Eqn. 
11071 Here, as in previous cases, we are free to structure the RG calculation for convenience. 
This easiest calculation is the second approach — in which only one constant of integration 
is actually renormalized — and our renormalized result is the same as our naive starting 
point. 

This simplified analysis (considering only the most secular terms) illustrates some of the 
pitfalls which can arise in applying RG to switchback problems. However, we must finish 
the O (e 2 ) analysis by considering all terms in Eqn. 11001 to understand the final nuances 
of this problem. There is a new complication when we attempt to renormalize all terms of 
Eqn. 11001 The final term, —B\e~ x ei(x), has the same kind of "splitting" ambiguity which 
we encountered in dealing with the doubly underlined term. 

We introduce our arbitrary position variable, r, which we want to choose so as to eliminate 
the secular term in x by replacing it with a divergence in r. In many cases, it is clear how 
to deal with the secular term. For example, a linear divergence — x — can be replaced 
with x — t + r. The final r will be absorbed into the renormalized constants of integration, 
and the x — r term (which is now considered non-secular), will ultimately disappear after 
renormalization. However the term —B\e^ x e\[x) is confusing. As seen above, there are two 
ways to "split" the B\e\(x)/2 term. There are four different ways to split e~ x ei(x). It may 
be replaced by any of the following: 

1. (e~ x — e~ T ) ei(x) + e~ T e\(x) 

2. e~ x e±(x) — e~ r ei(r) + e~ r ei(r) 

3. (e~ x — e~ T ) (ei(x) — e\{r)) + e~ T ei(x) + e~ x ei(x) — e _r ei(r) 

4. e~ x (ei(x) — ei(r)) + e~ x ei(r) 

All four of these options "cure" the divergent term (i.e. the secular term will vanish when 
we subsequently set r = x), and are equal to e~ x ei(x). If handled properly, any of these 
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options can lead to a valid renormalized solution. However, we will show that the fourth 
and final option is most natural, and results in the simplest algebra. 

How do we choose? The first consideration is subtle: The overall renormalized pertur- 
bative result must satisfy the governing equation (Eqn. l8U\) independently for each order in 
e. How we renormalize the O (e 1 ) divergences (Eqn. I9"5"j) has implications for O (e 2 ) calcu- 
lations. For example, in O (e 1 ) renormalization, there is an important difference between 
Eqn. [95] and Eqn. [851 The former has the additional term — eBi(r)ei(r). This term requires 
the presence of an additional O (e 2 ) term: e 2 e~ x B 2 (r)ei(r). Without this term the O (e 2 ) 
renormalized solution will not satisfy Eqn. [821 an d the renormalization procedure will yield 
an incorrect solution. We were able to gloss over this before because we were considering 
only the most secular term at O (e 2 ). 

Inspecting the four possible splittings enumerated above, we see that only the last two 
options provide the necessary e 2 e~ x Bf(r)ei(T) term, and can satisfy Eqn. [82] without con- 
trivances. 28 In examining both of these options, we split the e\{x) term for simplicity, 
as in the derivation of Eqn. 11071 29 Considering the third option first, our renormalized 
perturbation solution becomes: 

u{x) = l + e(B Q (r) + B 1 (r)(e 1 (x)-e 1 (r))) + e 2 ( K -B 1 (r)B (r)e- x - 
B?(r) (e- - e-) ( ei (x) - ei (r)) - ^(r) 2 {e\{x) - e 2 (r)) + 
2B 2 (r) ( ei {2x) - ei (2r)) ) + O (e 3 ) (108) 

As it must, this result satisfies Eqn. l8"2"1to O (e 2 ). By applying the RG condition (d T u(x) = 
0) to this equation, and grouping the resulting equation according to x dependence, we obtain 
a lengthy equation which can only be satisfied to O (e 3 ) \/x if: 

B[(ry = eB?(T) (109) 
e 2T TB' (r) = e 2r rei(r)fii(r) - e r J B 1 (r) - 3eB 2 (r) + e T eB 2 (r)ei(T) - e T reB^T^r) 
= e(eB 1 (r)+e T reBi(r)) 

28 The first two options can satisfy the governing equation if we carefully choose a different homogeneous 
solution at O (e 2 ) . With the proper non-zero choice of Co and C\ we can use the first two splittings 
enumerated, and they will result in an equivalent RG solution. 

29 In principle, each of the possible O (e 1 ) splittings could be paired with all possibilities at O (e 2 ), resulting 
in eight total possibilities. 
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Generally, no solution will exist, as we have two unknown functions and three differential 
equations. In this case, however, the first equation requires that: 

For this -Bi(t) solution, it is actually possible to satisfy the latter equations simultaneously 
to O (e 3 ): This occurs because the last equation is simply the lowest order of the second 
one. 30 There is another noteworthy point regarding the second part of Eqn. 11091 In all 
previous calculations, we discarded terms like e 2 -Bg(r), since B' (t) and B\{t) had to be of 
0(e 1 ). To solve these equations, however, B' (t) can not be Oie 1 ) (although B[(t) is). 
Solving for B , 

RW o r 2e + e°[3 1 + e° (2g - 1) eexfr) mn 
Bo(t)=(3 -J^ aie-erfr) ^ (1U) 

This solution, while valid, is cumbersome. Consider instead the fourth possible "split" 
enumerated above. Eqn. 11121 gives our renormalized perturbation solution, which satisfies 
Eqn. M 

u(x) = l + e(B (r) + B 1 (r)(e 1 (x)-e 1 (r))) + e 2 ^-B 1 (r)B (r)e- x - 
B\{T)e-* ( ei (x) - e x {r)) - l -B t {rf {e\{x) - e\{r)) + 
2Bl(r) (ei(2x) - e x {2r)) ) + O (e 3 ) (112) 



Applying the RG condition (d T u(x) = 0), and requiring that it be satisfied Vx, we obtain 
the following solutions for B (t), and -Bi(t): 

Bi(t) = /3i + C(e 3 ) (113a) 
B (r) = fa + PMr) + e f _#MM + 2/3 1 2 e 1 (2r)^) + O (e 3 ) (113b) 

Substituting these results into Eqn. 11121 and setting r = x, we obtain the final RG result, 
given by Eqn. 11141 

u(x) = l + e(j3 + /3 1 e 1 (x))+ (114) 
e 2 ( - Pifoe-* + 2f3 2 l e l {2x) - \(3le\{x) - P^e^x) ) + O (e 3 ) 



30 This can be seen explicitly by substituting Eqn. 11101 
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This is, of course, identical to our naive staring point, a happenstance we have seen several 
times previously. It is worth noting that the renormalized solutions obtained using Eqns. 
11101 and lllll are asymptotically equivalent to Eqn. 11141 

It may seem that we have needlessly digressed into the "terrible" problem. However, a 
clear-cut "best" strategy has emerged from our detailed exploration. Furthermore, we have 
identified — and resolved — a number of subtleties in the application of RG. Before applying 
these lessons to the problem of low Reynolds number flow past a cylinder, we summarize 
our conclusions. 

The "best" strategy is the one used to derive Eqn. I114[ a result which is identical to our 
naive solution (Eqn. HOOp . First, transform to the inner equation. Solve the O (e°) equation 
incompletely (obtaining just one constant of integration), which can then be set to satisfy 
the boundary condition at oo. This "trick" necessitates retention of integration constants 
at O (e 1 ), but results in computational simplifications (a non-linear RG equation) which are 
essential in dealing with the Navier-Stokes equations. 

At O (e 2 ), the homogeneous solution are identical to those at O (e 1 ). Consequently, the 
O (e 2 ) integration constants need not be retained, as we can meet the boundary conditions 
with the O (e 1 ) constants. We just pick a convenient particular solution. 

To apply RG to the terrible problem, we first "split" the secular terms. There are 
several ways to do this, even after requiring that the renormalized perturbation expansions 
satisfy the governing equations at each order. We can again choose for simplicity, bearing in 
mind that O (e 1 ) renormalization can impact O (e 2 ) calculations. It is easiest to apply the 
RG condition to the renormalized perturbation expansion, rather than applying it to the 
integration constants directly. In solving the resulting equation, we want solutions which 
are valid Vx. To solve the RG equation, care must be taken to satisfy several conditions 
simultaneously, and it cannot be assumed that our renormalized constants have a derivative 
of 0{e l ). 

Although there is quite a bit of flexibility in implementing the RG technique, our results 
are robust: Regardless of how we structure the calculation, our solutions agree to within 
an accuracy limited by the original naive perturbative solution; they are asymptotically 
equivalent. It is this robustness which makes RG a useful tool for the low Reynolds number 
problems, where the complexity of the Navier-Stokes equations will constrain our choices. 
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B. Flow past a cylinder 



1. Rescaling 



To solve Eqn. fusing RG techniques, we begin by transforming the problem to the Oseen 
variables. As in the terrible problem, to find a solution which is valid for all r, we need to 
analyze Eqn. [9] using a dominant balance argument. As it stands, different terms of Eqn. 
[9] will dominate in different regimes. 31 Looking for a rescaling of i/j and r which makes all 



terms of the same magnitude (more precise 



given in Eqn. 11151 ( Proudman and Pearson 



y, of the same order in R), yields the rescaling 



19571 ). 



p = Rr, 



* = Rip 



Transforming to these variables, Eqn. [9] becomes: 



The boundary conditions (Eqn. [TO]) become: 



P d(p,6) 



■*(p = R J 9) = J 



dp 



lim 



P=R 



p^oo p 



sin(8) 



;ii5) 



(116) 



117) 



2. Naive perturbation analysis 



The next step in obtaining the RG solution is to begin with the ansatz that the solution 
can be obtained from a perturbation expansion (Eqn. 11181) . 



*(p, e) = Mp, o) + R*i(j>, 0) + R 2 ^2(p, 0) + o (R 2 ) 



(118) 



31 i.e. the LHS, which is comprised of inertia! terms dominates for small |r| whereas at large |r| the viscous 
terms which comprise the RHS are of equal or greater importance. 
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Substituting Eqn. II 181 into Eqn. 11161 and collecting powers of R yields a series of equations 
which must be satisfied: 



o(R). v p Mp, o)--\ [-dffy- -^oe ) + \~WTp ~ "dim ) v ^ 



2 



<9# <9p <9p <96* 



V 2 *i 



5. solution 

The zeroth order part of Eqn. 11201 is the same as Eqn. 11161 and is equally hard to solve. 
But RG does not need a complete solution; we just need a starting point. We will begin with 
the equation which describes a uniform stream. This is analogous to the constant O (e°) 
solution in the "terrible" problem. 

A first integral to the O (R°) equation can be obtained by noting that any solutions of 
V 2 ,\I/o(P)$) = are also solutions of Eqn. 11201 This is Laplace's equation in cylindrical 
coordinates, and has the usual solution (assuming the potential is single- valued) : 

oo 

%(p, 6) = A + B \np + Y, {{A n p n + B nP - n ) sinnfl + (C n p n + D n p~ n ) cosnfl) (120) 

n=l 

We are only interested in solutions with the symmetry imposed by the uniform flow (Eqn. 
[TO]) . Hence Aq = Bq = C n = D n = 0. Furthermore, the boundary conditions at infinity 
require that A n = for n > 1. For simplicity at higher orders, we set C n = 0; this is not 
required, but these terms will simply re- appear at O (R 1 )- Finally set Ax — 1 to satisfy the 
boundary condition at oo (Eqn. II 17j) . As in the "terrible" problem, this is done for technical 
convenience, but will not change our results. We are left with the potential describing the 
uniform flow: 

* o (p,0)=psin(0) (121) 
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4- O (R 1 ) solution 



By substituting Eqn. 11211 into the O (R l ) governing equation, we obtain Eqn. 11221 



V^(P, 0) = cos(0)f - VJ*! (122) 



9 sin(0) 9 
' <9p p 89 

This equation is formally identical to Oseen's equation (Eqn. [371). albeit derived through a 
different argument. This is fortuitous, as its solutions are known (ITomotika and Aoil . Il950l ). 
Unfortunately, when working with stream functions, the solution can only be expressed as 
an infinite sum involving combinations of modified Bessel functions, K n , I n . 

The general solution can be obtained eithe r by following Tomotika or by using variation 
of parameters (IProudman and Pearson! . 119571 ). It is comprised of two parts, the first being 
a solution of Laplace's equation (as at O (R°))- The same considerations of symmetry and 
boundary conditions limit our solution: In Eqn. 11201 Aq = Bq = C n = D n = 0; A n = 0, if 
n > 1. Here, however, we retain the constants B n , and do not fix A±. This is analogous to 
what was done with the homogeneous terms at O (e 1 ) in the "terrible" problem. The second 
part of the general solution is analogous to a particular solution in the "terrible" problem, 
and can be obtained from Tomotika's solution (Eqn. |4*9|) . These two results are combined 
in Eqn. I123[ which will be the basis for our RG analysis. 

oo / oo \ 

*i(p,e) = A lP sm6 + B "P~ n + ^X mP ^ n {p/2) smn9 (123) 

n=l \ m=0 / 

Before discussing the application of RG to Eqn. 11231 it is worthwhile to discuss Eqn. 11221 
in general terms. Eqn. 11221 may be re-written as: 

^(V^cos^ + ^^V^O (124, 

We see explicitly that this equation is a linear operator (£) acting on ^i, and that the RHS 
is zero. This is the homogeneous Oseen equation. It is only because of our judicious choice 
of \l/o that we do not need to deal with the inhomogeneous counterpart, i.e. with a non-zero 
RHS. However, the inhomogeneous Oseen equation governs at all higher orders. This 
can be seen for O (R 2 ) from Eqn. 11201 

In general, the solutions to the inhomogeneous Oseen equation are found using the method 
of variation of parameters. It is worth exploring these solutions, as they provide some insight 
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into the structure of Eqn. HH We now solve Eqn. 1 1241 for a particular kind of inhomogeneity, 
one which can be written as a Fourier sine series. 32 We want to solve: 

oo 

£#i = ^ F n (p) sin n9 (125) 



n=l 



The substitution V 2 #i = e pcose/2 U 



'p, 6Q 33 , allows us to obtain th e first integral of Eqn. 



11251 This result is given by Eqn. 11261 (iProudman and Pearson! . Il957l ). 



(V 2 - U(p,e) = f^F n (p)sinn9 (126) 

Here F n (p) = e~ pcos9 ^ 2 F n (p). To solve for H(p, 9), begin by noting that the symmetry of 
the inhomogeneous terms implies that H(p, 9) can be written as a sine series. Consequently, 
substitute H(p, 9) = Y^=i9n(p) sinn9 into Eqn. 11261 to obtain: 

9n(p) + -9n(p) + 9n(p) = F n (p) (127) 



V4 ' P 2 / 

The fundamental solutions of this equation are K n (p/2), I n (p/2). Using variation of param- 
eters, the general solution of Eqn. 11271 may be written: 

9n{p) = -In (f) («n + jt\p)) + (0 (Pn + jt\p)) (128) 

Here, j} n \p) = J dppF n (p)K n (p/2), j} n \p) = J dppF n (p)I n (p/2), and a n , f3 n are con- 
stants. The next step is to undo our original transformation, and to solve the resulting 
equation: 



oo 

o , , p cos 8 v "\ , , 

V 2 ^i(p,#) = e^^„(p)sinn# (129) 

n=l 

oo 

= b n {p) sin n9 



n=l 



In this equation, b n (p) = J2m=i 9m{p) (I n -m (p/2) — I n +m (p/2)). We have the unfortunate 
happenstance that each b n depends on the all of the harmonics of the first integral. This is 
the origin of the nested sum (over m) in Tomotika's solution (Eqn. 1491) . 



32 The symmetry of the problem precludes the possibility of cosine terms in the governing equations for ty n , 
Vn > 1. 

33 V72 



V 2 4'i(p, 9) is the vorticity. 
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As before, symmetry will require that ^i(p, 6) be representable as a sine series: ^i(p, 6) = 
Ylm=i X m (p) sinm9. With this substitution we obtain (for each m), the radial component 
of Poisson's equation in cylindrical coordinates: 

X'i{p) + -X' m (p) - %X m {p) = b m (p) (130) 
p r z 

The fundamental solutions were discussed before in the context of Laplace's equation: p m , 
p~ m . As before, a particular integral is obtained through variation of parameters, and the 
general solution may be written: 

X n (p) = -P n (A n +l[ n \p)) + 1 (B n + ll n \p)) (131) 

Here z{ n) (p) = J dp - pb n {p) / \2np n ) , X^\p) = J dp - pb n {p)p n /{2n), and A n , B n are 
integration constants. 

It is useful to relate Eqn. 11311 to Tomotika's solution (Eqn. 149]) . There are four integration 
constants for each angular harmonic. Two are obvious: A n ,B n . The other two arise in the 
first integral (the vorticity solution), Eqn. 11281 However, every vorticity integration constant 
appears in each harmonic of Eqn. 11311 For example, one cannot uniquely assign ai and j3\ 
to the sin 9 harmonic of Eqn. 11311 However, if one considers n terms from Eqn. 11281 and n 
terms from Eqn. 11311 there will be 4ra integration constants — four per retained harmonic 
of Eqn. 11311 In passing we note that matched asymptotics workers avoid this problem by 
using the vorticity directly, and thereby simplify their treatment of boundary conditions. 
This approach does not work in conjunction with RG. 

It is mildly disconcerting to have four integration constants, as there are only three 
boundary conditions for each harmonic (Eqn. |117p . However, two of the constants — A n 
and a n — will be determined by the boundary conditions at infinity. This claim is not 
obvious, particularly since terms which are divergent prior to renormalization might not be 
present after the renormalization procedure. We outline here an argument which can be 
made rigorous. There are two kinds of divergences in Eqn. 11311 Terms which are secular as 
p — > 0, and terms which diverge too quickly as p — > oo. 34 

After renormalization, we will try to need to meet the boundary conditions (Eqn. 11171) . 
As in the case of the "terrible" problem, it will turn out that the simplest approach to 

34 To be precise, terms which diverge faster than p as p — > oo are problematic, and prevent satisfying the 
boundary conditions (Eqn. 11171) . 
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renormalization yields a renormalized perturbation solution which is the same as the naive 
series. Consider Eqn. 11311 The terms which are secular as p —>■ will not preclude satisfying 
the boundary conditions. Those which diverge too quickly as p — > oo, however, will conflict 
with Eqn. ITT71 

These terms must be eliminated by a suitable choice of integration constants. It turns 
out not to matter whether we do this before or after the renormalization procedure. For 
simplicity, we will do it before renormalizing. First, the coefficient of p n must vanish for all 
n > 1. This can happen, with an appropriate choice of A n , if 

lim l[ n) (p) -0(1) 

p— >oo 

For this requirement to be met, the coefficient of I n (p/2) in Eqn. 11281 must vanish (e.g., 
a n = Hindoo Ji(p)). It is always possible to choose a n appropriately, because the following 
condition is satisfied for all n: 

lim ^(p)- 0(1) 

p^OO 

In our problem this is true because F n (p) is based on solutions to the lower order governing 
equations. By construction, these are well-behaved as p —>■ oo. Therefore, for the inho- 
mogeneous Oseen equation under consideration (Eqn. 11261) . we see that two of the four 
integration constants — A n , a n — are needed to satisfy the boundary conditions at infinity. 

More specifically, the immediate problem requires us to consider the homogeneous Oseen's 
equation (Eqn. 11241) . and Tomotika's solution (Eqn. 09]). For this problem, F n (p) = 0, and 
the coefficient of I n (p/2) in Eqn. 11281 has no p dependence. So we simply choose a n such 
that this coefficient vanishes. Simplifying Eqn. 11281 we then have the following solution for 
the vorticity: 

oo 

V 2 *^) = e BS ^Y t K n (0 ((3 n )smn6 (132) 

n=l 

Since this solution for the vorticity is well-behaved as p —>■ oo, it follows that we can choose 
A n {n > 1) in Eqn. 11311 so that the coefficient of p n vanishes as p — > oo. We are left with 
the solution 

X n {p) = A n p5 n . l+ p n (l{ n \p) -Xi n) (oo)) +p~ n (B n +lt\p)) (133) 
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For the homogeneous Oseen's equation, l[ n \p) and (p) simplify to: 

4 n \ P ) = J dp^p- fcfMCr* (f) (in-m (f) - 4 +m (134) 

4 n) (p) = / dp^ P n (f> m i^ m (0 (/ n _ m (0 - I n+m (0) (135) 

This result is fundamentally the same as Tomotika's (Eqn. |49|) . However, his solution 
is more useful, as he accomplished the integrals in Eqn. 11341 What is the point of all 
this work? Firstly, the approach based on the variation of parameters may be applied to 
the inhomogeneous Oseen's equation, which must be solved for orders higher than O {R 1 )- 
Secondly, we see explicitly what happens to the two sets of integration constants a n and A n . 
Tomotika's solution has but two integration constants 35 — B n and (3 n . The other constants 
have already been chosen so as to satisfy the boundary conditions at oo. We have shown 
explicitly how they must be determined, and stated without proof that this may be done 
prior to renormalization. In short, we have explained why Eqn. 11231 is the appropriately 
general O {R 1 ) solution for our naive perturbation analysis. 

In addition to explaining why Tomotika's solution is a suitable starting point for 
RG, our analysis also connec ts with the O (R 1 ) solution of Proudman and Pearson 



( jProudman and Pearson! 119571 ) . We have shown that the vorticity must be well-behaved 
at p = oo if the overall solution is to satisfy the boundary conditions. 
a. Secular behavior 

Combining Eqns. 11211 11231 we begin the following naive solution: 

( p, 9) =psm(6) +R(Aipsm6 + jr ( B n p~ n + X m p$ m , n (0 ) smnOj + O (R 2 ) 

\ n=l \ m=0 / / 

(136) 

Although intimidating, this is conceptually equivalent to Eqn. [85] (in the terrible problem). 
The first step in our analysis is identifying which terms are divergent. As explained above, 
Eqn. 11361 is specifically constructed to be of O (p 1 ) as p — > oo. In fact, only the O (R°) and 
Ai terms matter at large p. As p —>■ 0, however, many other terms in Eqn. 11361 diverge. All 
of the B n terms diverge. Most of the $ mjn (p) terms are also secular. 



:>,-> 



There is also A\, but that is a special case. 
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Rather than enumerating and sorting through the different divergences, we simply treat 
the problem abstractly. Eqn. 11361 can be rewritten as: 



^ { p, 9)=p sm(9) + R (A lP sin 9 + K(p, 9; {B,}; {Xj}) + S(p, 9; {B m }; {X n })) (137) 
Here, S includes the terms which are secular as p — > 0, and 1Z includes regular terms. 

5. Renormalization 

Equation 11371 is renormalized just like the terrible problem. We begin with the renormal- 
ized perturbation expansion given in Eqn. 11381 Note that we are not specifying the details 
of which terms are secular, or how we are "splitting" these terms. The only term we are 
explicitly considering is A\. This is a trick built on consideration of the terrible problem. 
Our "best" solution (Eqn. II 141) to that problem was built on the renormalization of just 
one constant, B Q in Eqn. I113al Essentially, we will repeat that procedure here, using A\ as 
that constant. 

9) =psin(9) + R^A l (r)psm9 + K(p, 9; {B t (r)}; {Xj(r)}) + (138) 
S(p, 9- {B m (r)}; {X n (r)}) - S(r, 9; {B m (r)}; {X n (r)}) + O (R 2 ) 



We will now apply the RG condition — d T ^(p, 9) = — to this equation. Accomplishing 
this in complete generality is difficult. However, using our experience from the terrible 
problem, we can see that this is not necessary. The RG condition may be satisfied as 
follows: First, suppose that X n (r) = O (R 2 ) Vn, B m {f) = O (R 2 ) Vm. These equations 
are satisfied by X n (r) = Xm B m {r) = /3 m . Substituting these results into Eqn. I138j and 
applying the RG condition results in: 

= R (A' 1 (T)psm9- S'(r,9; {f3 m }; { Xn })) (139) 

This is easily solved for A\{r). 

S(T,9;{M;{xn}) +ai (14Q) 

psin^ 

We have explicitly validated our supposition that {X n {j)} and {B m (r)} can be constants. 
With this supposition, we have shown that the RG condition applied to Eqn. 11381 can be 
satisfied with an appropriate choice of A\{r). We have satisfied the RG condition through 
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clever tricks derived from our experience with the terrible problem. However, this solution 
is entirely valid, and our experience with the terrible problem has shown us that more 
complicated solutions are asymptotically equivalent. 

Substituting Eqn. I14UI into Eqn. 1138} and setting r = p, we obtain our renormalized 
solution: 

^ {P , 9)=p sm(9) + R (a lP sin 6 + TZ(p, 6; {#}; { Xj }) + S(p, 6; {(3 m }; { Xn })) (141) 

By now it should not be surprising that this is the same equation as our naive pertur- 
bation solution (Eqn. [1371) . and by extension the same solution obtained by Tomotika 



( ITomotika and Aoil . Il950l ). As in the case of the terrible problem, however, we now know 
that this is a uniformly valid approximation. We now may choose the integration constants 
to satisfy the boundary conditions, and then calculate the drag coefficient. 
a. Truncation 

Unfortunately, there are infinitely many integration constants, and it is impossible to 
apply the boundary conditions to our renormalized solution (or Eqn. 11361) . To progress 
further, we must make the sam e sort of uncontrolled approxim ations made by previous 



workers (IProudman and Pearson! . 119571 ; ITomotika and Aoil . Il950f ) . 
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Our approximation consists of a careful truncation, in both m and n, of the series in Eqn. 
11361 There are two important points to consider. First is the sin 9 symmetry of the overall 
problem: terms proportional to sm9 reflect the symmetries exhibited by the uniform flow 
which are imposed on our solution via the boundary conditions at infinity. The importance 
of this harmonic is further seen in Eqn. [2TJ Only the coefficient of sin9 will be needed for 
the computation of Co- 

Secondly we recall that the remaining boundary conditions are imposed at the surface 
of the sphere, at p = R in Oseen coordinates. When applying the boundary conditions, 
terms which are secular as p — > will therefore be most important. Specifically, we cannot 
truncate any terms which are divergent, although we are at liberty to set their coefficients 
equal to zero. 

These considerations allow exactly one solution. First, set all B n = n > 1. Secondly, 



36 Kaplun was able to avoid this difficulty by using the velocity field instead of stream functions, although 
his approach brings other problems: the solution cannot be expressed in closed form, and must be ap- 



proximated to apply the boundary conditions (see section [lI.B.3p . 
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n = 


1 2 3 4 


*£Ld(p.*) 


O (i? 3 log #) O {R 2 ) (R 1 ) O (R°) 
0(R 2 \ogR) 0(R 1 ) 0(R 1 ) 0(R- 1 ) 



TABLE IV Relative importance of discarded terms at p = R. 



set all X m = m > 0. We retain three coefficients: Ai,Bx,X , which will permit the 
boundary conditions to be satisfied for the sin# harmonic. What about the higher harmon- 
ics? These terms are truncated in an uncontrolled approximation. However, as we will show, 
the discarded terms are O (R 3 log R) or higher at the surface of the sphere. They are regular 
terms, and thus negligible in comparison to the secular terms retained (which are O (R^ 1 )). 

Now, suppose we follow Tomotika, and try to extend this approach, by retaining a few 
more terms. The next step would be to retain the B 2 ,Xi terms, and to try to satisfy the 
boundary conditions for the sin 28 harmonic. As before, all the higher B n ,X m are set to 
zero. Why not include the next harmonic or two? 

The answer lies in the terms we discard. If we satisfy the boundary conditions at p = R 
for the first n harmonics, we must retain the coefficients Xq, ...,X n — 1. To minimize the 
amount of truncation we do, first set X m = for Vm > n — 1 and Bk = for VA; > n. What, 
then, is the form of the terms which are discarded from our solution? 

too n—1 \ 
^2XA, k (p/2)p S inke\ (142) 
k=n+l m=0 / 

^discard (/)' 9) is largest as p — > 0, and will be most important at p = R, on the surface of 
the cylinder. If we retain only the n — 1 harmonic, ^dis car( j (p, 9) ~ 0(R 3 \ogR). Since 
we are only working to O {R 1 ), this is fine. We must also consider the derivative, since we 
want to satisfy all of the boundary conditions (Eqn. [101) to the same order. ^ diS c ar d(P! ~ 
O (R 2 logR) Therefore, in the case where we retain only the sin# harmonic, the discarded 
terms are negligible, as we are working to 0(R 1 ). 37 When we retain higher harmonics, 
everything changes. Table IIVI shows the magnitude of the discarded terms at p = R for the 
first four harmonics. 



37 This argument is somewhat simplistic: The neglected terms also contribute, when meeting the boundary 
conditions, to the values of the retained coefficients, i.e. All non-zero X m affect Xq. But these are lower 
order effects. 
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From Table IIV} we see immediately that to retain sin 29 harmonics, we must have an 
error in our derivative boundary condition of O (R 1 ) — the order to which we are trying 
to work. If we retain higher harmonics, this situation gets worse. First we have an O {R 1 ) 
error in the stream function itself, and then we begin to have errors which are divergent in 
R\ For n > 4, both ^^l cavd {p, 9) and ^^ C3X ^{p-i 9) are increasingly divergent functions of R. 

Since it is in practice impossible to fit the boundary conditions to Eqn. 11361 we must 
truncate the series expansion. We have shown that there is only one truncation consistent 
with both the symmetry requirements of the problem and the demand that we satisfy the 
boundary conditions to O (R 1 )- 

9)=p sm(0) + R (a iP + B IP - 1 + Xop^o,! (0 ) sin 9 (143) 

This result is identical to Proudman's O {R 1 ) result for the Oseen stream function 



( iProudman and Pearsonl . 119571 ). However, he arrives at this result by considering matching 



requirements with the O (R°) Stokes expansion and by imposing sin 9 symmetry on the first 
integral (Eqn. I130p . Our approach arrives at the same conclusion, but without the need 
for asymptotic matching or the two expansions it requires. Moreover, we did not need the 
expertise and finesse which matched asymptotics workers needed to deduce the unusual form 
of their expansions (e.g., the 1/logi? ter m in Eqn. 155]) . Finally, w e note that Tomotika's 
numerical results support our truncation (ITomotika and Aoil . Il950l ). 
b. Meeting boundary conditions 

It is straightforward to apply the boundary conditions (Eqn. [10]) to Eqn. 11431 To satisfy 
the condition at infinity, A\ — 0. The other two requirements are met by the following 
choice of coefficients: 

-R 2 $'{R/2) t s 

Bi = , ; \ ' - , i \ 144 

-4 

X ° = R(41> 0il (R/2) + R& 0il (R/2)) {U5) 

Notice that we are using the Oseen stream function. The Stokes' function is related by: 
ip(r,9) = ty(rR,9)/R. Putting everything together, we have the new result given by Eqn. 



[T461 



ty { p,9) = psm(9) + R 



-R^' 0>1 (R/2) 




(146) 



4®o,i(R/2) + R& 0>1 (R/2) 



-4 



p$o,i(p/2) sintf 



i?(4$ , 1 (i?/2) + i?$o il (i?/2)) 



Remember that although our truncated solution satisfies the boundary conditions exactly, 
it only satisfies the governing equations approximately. 

6. Calculating the drag coefficient 

We now transform Eqn. 11461 into Stokes' coordinates, and substitute the result into Eqn. 
I2T1 38 We thereby obtain a new result for Co, given by Eqn. 11471 



This result is plotted in Figure [T4"l where it is compared against the principal results of Oseen 
theory, matched asymptotic theory, and experiments. When compared asymptotically, all 
of these theoretical predictions agree. At small but not infinitesimal Reynolds number, the 
largest difference is seen between Kaplun's second order result and the first order predictions, 
including Eqn. 11471 As explained previously, current experimental data cannot determine 
whether Kaplun's second order matched asymptotics solution is actually superior. 

The RG result lies among the first order predictions. Fundamentally, the RG calculation 
begins with an equation similar to Oseen's, so this is not too surprising. Within this group 
Eqn. 11471 performs very well, and is only slightly bettered by Imai's prediction (Eqn. I5T1) . 
These two results are very close over the range < R < 1. 

The real strength of Eqn. 11471 can be seen in in Figure [151 As the Reynolds number 
increases beyond R — 1 , all other theories begin to behave pathologically. They diverge from 
experimental measurements and behave non-physically (e.g., a negative drag coefficient). 
The RG prediction suffers from none of these problems; it is well behaved for all R. As it is 
still based on a perturbative solution, it does become less accurate as R increases. 

38 Or, alternatively, into Eqns. [HHH and 1201 
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TV 



(-12^(^/2) + R (6*Q tl (R/2) + R^jR/2))) 
8$ 0tl (R/2) + 2R& 0tl (R/2) 



(147) 
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FIG. 14 (Color online) Drag on cylinder, comparing RG predictions to other theories at low R. 



C. Flow past a sphere 



1. Reseating 

Our analysis of low Reynolds number flow past a sphere closely follows both the cylinder 
problem and the terrible problem. We omit redundant explanations. As before, the first step 
is a rescaling of both r and ip — the transformation into Oseen coordinates. A dominant 
balance analysis identifies the rescaling given in Eqn. 11481 

p = Rr, * = R 2 ip (148) 

In Oseen variables, the governing equation (Eqn. [1~2|) becomes: 

1 fd(ty(p,u),D*-9(p,u)) „ \ 
D$*(p, fi) = j 2 ( d[p I) + 2D >^ ri L >*(P, A*) J (149) 



where 

= pnsfl D 2 = r . 

dp 2 p 2 dp 2 ' 1 — p 2 dp p dp 

90 



, scose , fl. s * + iZ^|L, L , s _JL» + i» (150) 




The boundary conditions (Eqn. [T^I) transform into: 



*(p = J R,/i) = 0, 



<9p 



0. 



p=R 



lim 



(151) 
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2. Naive perturbation analysis 



We continue by substituting our naive perturbation assumption (Eqn. 11181) into Eqn. 
11491 and then collecting powers of R. 

O (If) : D tM P, p) = ^^(P,g,^0(p^)) + 2D , Mp) ^ (152a) 
<>(*): = - ^ ^ + + 

2 (£)J* ^*i + ^i^o) j (152b) 
,/ p21 n4xW v 1 ( d(* ,Dl* 2 ) d^D^} d(* 2 ,D* p * ) 

2 (Dj* ^*2 + D^L^ + £?* 2 L p *o) J ( 152c ) 



5. solution 

As seen with both the cylinder problem and the terrible problem, Eqn. I152al is the same as 
the original governing equation (Eqn. |149p . As before, we proceed using an incomplete solu- 
tion for \l/o : the uniform stream which describes flow far from any disturbances. Analogously 
to the cylinder, we notice that Eqn. I152al is satisfied if ^o(p, p) obeys D^ (p,p) = 0. The 
general solution of this equation which also satisfies the appropriate symmetry requirement 
(* (p,/x = ±1) = 0) is given by Eqn. QSSJ 

oo 

Mp, p) = J2 ( A -p n+1 + B ^ n ) (153) 

n=0 

Here Q n {p) is defined as in Eqn. SSI Following the analysis used for the cylinder, we set 
all of the coefficients to zero, excepting A\ = —1/2. This choice of A\ satisfies the uniform 
stream boundary condition (Eqn. 11511) at p = oo. We thereby obtain: 

$ (m) = 7QiM (154) 
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4- O (R 1 ) solution 



Substituting Eqn. 11541 into Eqn. I152b| we obtain Eqn. 11551 

DtMp, a*) = ( + v-f) d 2 p Mp, A*) 



(155) 



This result is also derived in matched asymptotic analysis, and is formally identical to the 
Oseen equation for a sphere (Eqn. |3"51) . Stru cturally, this pro blem is similar to what we have 



seen previously, and is solved in two steps ([Goldstein 
DH> X = e^/ 2 $(p,p) to obtain Eqn. [T561 39 



19291 ) . First use the transformation 



L> p 2 --)$(p,p) = 



This may be solved to obtain the first integral: 



71=1 



P 

•5 V2 



B r , 



p\2 



Qn{p) 



(156) 



(157) 



As in the case of the cylinder, the inhomogeneous term on the RHS of Eqn. 11571 consists 
of integration constants which multiply the two modified Bessel functions. We are beset by 
the same considerations, which (properly speaking) must be resolved by applying boundary 
conditions (Eqn. 1 1 51 j) to the renormalized solution. Following the same arguments given for 
the cylinder, we set the coefficients B n = 0, which will later make it possible to satisfy the 
boundary conditions at infinity. 

Comple t ing t he second integration is difficult, but was accomplished by Goldstein 
( iGoldsteinl . Il929l ). The requisite solution is essentially the second term in Eqn. H6j 



^\p,6) = A 1 p 2 Q 1 (p) + B nP~ n + ^X mP 2 $ m , n (p/2) Q n {p) (158) 

n=l \ m=0 / 

Note that we have omitted the terms A n r n Q n (p) which diverge too quickly at infinity (this 
was also done for the cylinder). 

Alternatively, one may simplify the series in Eqn. 1157} by retaining only the n — 1 term 
(setting all other A n — 0). It is then possible to complete the second integration with a 
closed form solution: 



v|>S b) (p, 6) = A 1 p 2 Q 1 (p) +Ai(l + p) (l - e -ioM) + ^ B n p- n Q n {p) 



(159) 



71=1 



:j>) 



Z3p*i(p, n) is the vorticity. 
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As before, we neglect the A n r n Q n (p) solutions. This is essentially Oseen's solution (Eqn. 
[38]) . expressed in the appropriate variables and with undetermined coefficients. 

We therefore have two solutions (Eqns. 11581 1 1591) which can be used for vE^. For the 
moment, we will consider both. We will later demonstrate that the former is the preferred 
choice by considering boundary conditions. 

5. Secular behavior 

We consider our O {R l ) naive solution abstractly: 

^(p,/i) = -p 2 Q 1 {p) + R i^A lP 2 QM + f^B n p~ n Q n {p) + • • - j + O (R 2 ) (160) 

This generic form encompasses both Eqn. 11591 and Eqn. 11581 It also possesses two key 
similarities with both the terrible and the cylinder problems. First, there is a term at O (R 1 ) 
which is a multiple of the O (R°) solution (A 1 p 2 Q 1 (p)). Secondly, the secular behavior in 
our naive solution occurs at the same order as the integration constants which we hope to 
renormalize. 40 This fact is in essence related to equations like Eqn. [891 which must be solved 
iteratively. We avoided that kind of RG equation by introducing the constant which could 
have been associated with the O (R°) solution at 0(R l ). But renormalizing divergences 
into integration constants at the same order limits the ability of RG to "re-sum" our naive 
series. In all of these cases, the real power of RG techniques could be seen by extending our 
analysis to O (R 2 ). 

Because of the similarities between Eqn. 11601 and Eqn. 11361 we can tackle this problem 
in a manner formally the same as the cylinder. By construction, Eqn. 11601 is O (p 2 ) as 
p — > oo. Hence the only terms with problematic secular behavior occurs in the limit p — > 0. 
As before, these divergences need not even be explicitly identified. We write: 

( p, y) = -p 2 Qi{p) + R (A 1 p 2 Q 1 (p) + K{p, p; {B,}- {Xj}) + S{p, p; {B m }; {X n })) (161) 

Here, S includes the terms which are secular as p — >• 0, and 1Z includes regular terms. 

40 These secular terms are not written explicitly in Eqn. 11601 They can be found in Eqns. 11591 and 11581 



94 



6. Renormalization 



Eqn. 11611 is only cosmetically different from Eqn. 11371 Renormalizing the two equa- 
tions can proceed in exactly the same fashion. Therefore, we may immediately write the 
renormalized solution: 

y) = -p 2 Qi(p) + R («ip 2 Qi(/i) + K(p, 6; {&}; { Xj }) + S(p, 6; {(3 m }; { Xn })) (162) 

This is, of course, the same solution from which we began. As in the previous two problems, 
we now know that it is a uniformly valid solution, and turn to the application of the boundary 
conditions. 

7. Meeting the boundary conditions 

We have two possible solutions for ^i(p, fi). Considering the boundary conditions on 
the surface of the sphere (Eqn. 11511) will demonstrate why Eqn. 11581 is preferential. Eqn. 
11591 can never satisfy the two requirements for all of the angular harmonics. Expanding the 
exponential term, we see that although it has but one integration constant, it contributes to 
all of the powers of p. The second solution, Eqn. 11581 can meet both of the boundary con- 
ditions — in principle. However, as in the case of the cylinder, this is practically impossible, 
and we must consider truncating our solution. 

It is clear that we will need to approximate our solutions in order to apply the boundary 
conditions. Our procedure is governed by the following considerations. First, we demand 
that our approximate solution satisfy the boundary conditions as accurately as possible. 
This requirement is necessary because our goal is to calculate the drag coefficient, Cr>, a 
calculation which is done by evaluating quantities derived from the stream function on the 
surface of the sphere. Hence it is necessary that the stream function be as accurate as 
possible in that regime. Secondly, we want the difference between our modified solution and 
the exact solution (one which satisfies the governing equations) to be as small as possible. 

a. Oseen's solution 

First, consider trying to satisfy these requirements starting from Eqn. 11591 Although 
this is the less general solution to Oseen's equation, we consider Oseen's solution because 
of (1) its historical importance, including widespread use as a starting point for matched 
asymptotics work and (2) the appealling simplicity of a closed-form solution. 
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We combine Eqns. 11591 and 11541 to begin from the solution: ty(p,p) = fy (p,p) + 
Ri&i (p, p). Since we are interested in the solution near the surface of the sphere (p = R), 
and because there is no other way to determine the integration constants, we expand the 
exponential in that vicinity. Retaining terms up to O {Rp 1 ) ~ O (p 2 ), we obtain: 



y) = (V + R {A lP 2 - A lP )) QM + R^ B nP ~ n Q n (p) 



(163) 



?1=1 



The boundary conditions are satisfied if B n — Vn > 1, A\ — 0, Ai — —3/2, and B\ = 
—R 2 /2. In passin g, we note tha t substituting these values into Eqn. 11591 reproduces Oseen's 



original solution (jOseen 
obtaining: 



1910) . Continuing, we substitute these values into Eqn. 11631 

R 3 ' 



V + 3Rp 



(164) 



2 2p 

This is nothing more than Stokes' solution (Eqn. 128]) . albeit expressed in Oseen variables. 
Consequently, when substituted into Eqns. HHEU and [251 Eqn. 11641 reproduces Co = 6-7T / R. 

How accurate is our approximate solution? The difference between Eqn. 11641 and Eqn. 
11591 is given by: 

A\I> = --R (1 + p) (-2 + 2e-3^ 1 -") + p (1 - p)^j (165) 

At the surface of the sphere (p = R), this equates to an O {R 3 ) error in the stream function, 
and an O (R 2 ) error in the derivative. That is entirely acceptable. However, at large p, A^> 
grows unbounded, being of O (p 1 ). This is the fundamental problem with the solution given 
by Eqn. 11641 By beginning from Eqn. 11581 we can avoid this difficulty 

It is at first a little dis concert i ng th at Oseen used his solution to obtain the next approxi- 
mation to Cd (Eqn. [39]) (jOseenl . Il913l ). How can our results be worse? As noted previously, 
"Strictly, Oseen's method gives only the leading term ... and is s carcely to be counted as 



super ior to Stokes' method for the purpose of obtaining the drag." (IProudman and Pearson 



19571 ) 



b. Goldstein's solution 

We now apply the boundary conditions to Eqn. 11581 By starting from the more general 
solution to Oseen's equation, we can remedy the difficulties encountered above. This analysis 
will be very similar to the truncation performed on Tomotika's solution for the cylinder 
problem. 

We combine Eqns. 11581 and 11541 to begin from the solution: *ff(p,p) = *ff (p,p,) + 
(p, p). As with the cylinder, we will approximate the full solution by truncating the 
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n = 


12 3 4 


*c£card(^) 
^discard (P' p) 


O (i? 3 ) O {R 2 ) O (R 1 ) O (R°) 
O (R 2 ) O (R 1 ) O (R 1 ) O (R- 1 ) 



TABLE V Importance of discarded terms at p = R. 

series in both m and n. Our first consideration is again symmetry: The uniform flow imposes 
a sin#, or Qi(fi) symmetry on the problem. Hence we must retain the n = 1 term in Eqn. 
11581 The importance of this term is clearly seen from Eqn. [26l Only the coefficient of Q\ (/x) 
is needed to calculate the drag if the stream function satisfies the boundary conditions. 

As in the case of the cylinder, if we retain n harmonics, we must retain m = n — 1 
terms in the second sum (the sum over m) in order to meet both boundary conditions. 
To minimize the error introduced by our approximations we set all other B n , X m equal to 
zero. The remaining terms, those which would violate the boundary conditions and must 
be truncated, are then given by Eqn. 11661 

/ oo n— 1 \ 

*2ca r d(P> P) = R[ E E X m ^ m , k {p/2)p 2 QM (166) 

\k=n+l m=0 / 

We want to estimate the magnitude of the error in our approximation, both overall and 
at the surface (the error in the boundary conditions). The error is given by Eqn. 11661 First, 
we calculate the magnitude of both ^^scard(P' P) anc ^ ^ s derivative at the surface (p = R) 
with n retained harmonics. The results are given in Table PVT . 

From Table El we see that to retain the (^(a*) harmonics, we must have an error in our 
derivative boundary condition of O {R 1 ) — the order to which we are trying to work. If we 
retain higher harmonics, this situation gets worse. 

Since it is in practice impossible to fit the boundary conditions to all harmonics, we must 
truncate the series expansion. We see that there is only one truncation consistent with both 
the symmetry requirements of the problem and the demand that we satisfy the boundary 
conditions to O (R 1 ): 

V(p, p) = VQi(aO + R (A lP 2 + B x p- X + X $o,i(p/2)p 2 ) + O (R 2 ) (167) 

We also must consider the overall error, e.g., how big can ^iLird 0°> P) Although, at 
the surface of the sphere, Eqn. 11671 is no better than Eqn. I164[ it is superior for p ^ R. 
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The magnitude of the error is maximized as p — > oo. It can be shown by Taylor expansion 
(separately accounting for m = 0,m = n, etc.) that $ m , n (^ — > oo) ~ x~ 2 . Therefore, 



lim *2Ld(p,A*)=0(i2 1 ) 



Although this is somewhat unsatisfactory, this solution does not suffer from the same short- 
comings as Eqn. 11631 The error remains bounded. 

Eqn. 11671 will satisfy the boundary conditions (Eqn. I15ip if A\ = and 

X ° = 6R$ 0tl (R/2) + ^$0,1(^/2) (168) 

1 6f , 1 (iJ/2) + M;, 1 (iJ/2) 1 J 

As in the case of the cylinder, the resulting stream function satisfies the boundary conditions 
exactly, and the governing equations approximately. Our final solution is: 

/ i? 3 $n 1 (R/2) 

«to /.) = + «( Moil(fl/2) %;;^ /2) ^' + a™) 

For reference, 



3tt / 2 1 e~ 2a: 
$0,1 (x) = --^ 2- - + — - 



8. Calculating the drag coefficient 

We calculated the drag coefficient by substituting Eqn. II 701 into Eqn. [261 giving this new 
rcsu.lt ■ 

c = vr (-16j M (j/2) + g (8<&; tl (fl/2) + ggi(j/2))) 
D 2(6$ 0|1 (Ji/2) + fl$; il (iZ/2)) 

This can be expressed in terms of more conventional functions by substituting for <3> 0) i(x), 

resulting in the drag coefficient given by Eqn. 11721 

4tt (24 + 24i? + 8R 2 + R 3 + Ae R (R 2 - 6)) 
° D = R(2(R + l) + e R (R 2 -2)) (172) 

This result is plotted in Figure [161 where it is compared against the principal results of 
Oseen theory, matched asymptotic theory, numerical results, and experiments. As R — > 0, 
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on a sphere, comparing RG to other theories (jDennis and Walkerl . 



1970; 



Maxworthy 



1965|). 



there is excellent agreement. At small but non-infinitesimal Reynolds numbers, RG is nearly 
identical to Oseen's prediction (Eqn. [39]), which is disappointing. It is surprising that 
Goldstein's result is better than the RG result, as they are calculations of the same order in 
R, and are a series approximation. That the matched asymptotics predictions are superior 
is not surprising; Chester and Breach's result began with a much higher order perturbative 
approximation. If a higher order RG calculation were possible, RG ought to be better than 
the same order matched asymptotics prediction. 

As in the case of the cylinder, the real strength of Eqn. 11721 can be seen as the Reynolds 
number increases. Figure [T7| demonstrates that all other theories diverge from experimental 
measurements for R > 1. This is an unavoidable aspect of their structure and derivation - 
they are only valid asymptotically. The RG prediction suffers from none of these problems. 
Eqn. 11721 is well behaved for all R, although it does become less accurate at larger Reynolds 
numbers. 
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IV. CONCLUSIONS 

We have devoted a substantial effort to the historical problem of calculating the drag 
coefficient for flow around a cylinder and a sphere at low Reynolds number. We report 
four principal accomplishments. First, we have untangled over 150 years of diffuse, confus- 
ing, and sometimes contradictory experimental, numerical, and theoretical results. We have 
expressed all important previous work within a consistent mathematical framework, and 
explained the approximations and assumptions which have gone into previous calculations. 
Moreover, by plotting experimental results and theoretical predictions with the leading order 
divergence removed (an idea originally due to Maxworthy), we have consistently and criti- 
cally compared all available measurements. There are no other such exhaustive comparative 
reviews available in the existing literature. 

Secondly, we have extended traditional matched asymptotics calculations. We advance 



100 



and justify the idea that uniformly valid approximations, not the Stokes or Oseen expansions, 
should be used to calculate derivative quantities such as Cd- By combining this approach 
with previously published matched asymptotics results, we obtain new results for the drag 
coefficients. These results systematically improve on published drag coefficients, which relied 
only on the Stokes expansion. This methodology also resolved a problem in the existing 
liter ature: the mos t accurate calculations for a cylinder, due to Skinner, had failed to improve 



Cd (ISkinnerl . 



19751 ). When treated via a uniformly valid approximation, our new result based 



on Skinner's solutions betters all matched asymptotics predictions. 

We have also explored the structure and subtleties involved in applying re normalization 



group techniques to the "terr ible" problem posed by Hinch and Lagerstrom (IHinchl. 



1991 



Lagerstrom and Castenl . Il972l ) . This problem, previously solved by Chen et al. fjChen et al. 



19961 ) . contains a rich and henceforth unexplored collection of hidden subtleties. We exhaus- 



tively examined all possible complications which can arise while solving this problem with 
the renormalization group. To treat some of these possibilities, we identified and imple- 
mented a new constraint on the RG calculation; the renormalized perturbation solution 
itself, not just the expansion on which it is based, must satisfy the governing equations to 
the appropriate order in e. While this had been done implicitly in previous calculations, we 
had to deal with it explicitly (e.g., by appropriate choices of homogeneous solutions). In 
the process of doing so, we obtained several new second order approximate solutions to the 
"terrible" problem, and demonstrated their equivalence. 

The work with the "terrible" problem laid the foundation for our most significant new 
calculation. In close analogy with the "terrible" problem, we used the RG to derive new 
results for the drag coefficients for both a sphere and a cylinder (Eqns. 11721 and 11471 
respectively). These new results agree asymptotically with previous theoretical predictions, 
but greatly surpass them at larger R. Other theories diverge pathologically, while the results 
from the RG calculation remain well behaved. 

We demonstrated that these new techniques could reproduce and improve upon the re- 
sults of matched asymptotics — when applied to the very problem which that discipline was 
created to solve! Matched asymptotics requires the use of two ingenious and intricate expan- 
sions, replete with strange terms (like RlogR) which must be introduced while solving the 
problem via a painful iterative process. RG requires only a single generic expansion, which 
can always be written down a priori, even in complicated singular perturbation problems 
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with boundary layers. It therefore gives rise to a much more economical solution, requiring 
half the work and yielding a superior result. It is hoped that demonstrating of the utility 
of these techniques on this historical problem will result in increased interest and further 
application of renormalization group techniques in fluid mechanics. 
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